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XV. On a Class of Differential Equations, including those which occur in Dyna- 
mical Problems.— Vart II. By W. F. Donkin, M.J., F.R.S., F.R.A.S., Savilian 
Prof essor of Astronomy in the University of Oxford. 

Received February 17, — Read March 22, 1855. 

The following paper forms the continuation and conclusion of one on the same 
subject presented to the Royal Society last year, and printed in the Philosophical 
Transactions for 1854. I have however put it, as far as possible, in such a form as 
to be independently intelligible. 

The fourth Section (the first of this Part) contains a recapitulation of some of the 
most important results of the former Part, in the form of seven theorems, here enun- 
ciated without demonstration. 

In the fifth Section the method of the variation of elements is treated under that 
aspect which belongs to it in connexion with the general subject. It is applied, by 
way of example, to deduce the expressions for the variations of the elliptic elements 
of a planet's orbit from the results of art. 30 (Part I.), on undisturbed elliptic mo- 
tion; this example was chosen, partly because the resulting expressions are required 
in a future section, and partly for the sake of incidentally calling attention to a fal- 
lacy which has been, perhaps, often committed, and certainly seldom noticed. The 
same method, under a slightly different and possibly new point of view, is applied, as 
a second example, to the problem of the motion of a free simple pendulum, omitting 
the effect of the earth's rotation. I believe the methods of this paper might be advan- 
tageously employed in the treatment of that general form of the problem of a free 
pendulum which has been considered by Professor Hansen in his Prize Essay. I was 
unwilling, however, to attempt what might have turned out to be merely an uncon- 
scious plagiarism, without having seen the Essay in question, of which I only suc- 
ceeded in obtaining a copy on the day of writing this preface. As I now perceive 
that the investigation would be quite independent, I hope to enter upon it at some 
future time. 

The sixth Section contains some general theorems concerning the transformation 
of systems of differential equations of the form considered in this paper, by the sub- 
stitution of new variables. The most important case consists in the transformation 
from fixed to moving axes of coordinates, in dynamical problems. Some of the 
results are, I think, interesting, and perhaps new. 

The seventh and last Section contains an application of the preceding theorems, in 
connexion with the variation of elements, to the transformation of the differential 
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equations of the planetary theory. This investigation, if interesting at all, will pro- 
bably be so to the mathematician rather than to the astronomer. I think, however, 
that if the theories of physical astronomy were more frequently treated rigorously 
and symmetrically, apart from any approximate integrations ; and if, when the latter 
are introduced, more care were taken to give a clear and exact view of the nature of 
the reasoning employed, it might be possible to draw the attention and secure the 
cooperation of a class of mathematicians who now may well be excused, if, after a 
slight trial, they turn from the subject in disgust, and prefer to expatiate in those 
beautiful fields of speculation which are offered to them by other branches of modern 
geometry and analysis. 

The contents of the two last Sections are more or less closely connected with the 
subjects of various memoirs by other writers, especially Professor Hansen and the 
Rev. B. Bronwin. I cannot pretend to that degree of acquaintance with them which 
would enable me to give an exact statement of the amount of novelty to be found 
in my own researches. I believe it is enough to justify me in offering them to the 
Society ; beyond this I make no claim. 

Oxford, Feb. 15, 1855. 

Section IV. 

49. The following theorems were demonstrated in the former part of this essay, 
and are I'ecapitulated here for convenience of reference. (As before, total differen- 
tiation with respect to the independent variable t will, in general, be denoted by 
accents, which will be usedybr no other purpose.) 

Theorem I. — If X be a function of w variables Xj, .x'a, ...«„, and if 3/1,3/2? •••y« be w 
other variables connected with the former by the n equations 

dX_ dK_ dX_ 

then will the values oi x^, x^, ....*„, expressed by means of these equations in terms 

of «!,.... w„, be of the form 

_dY _dY _dY 

and Up be any other quantity explicitly contained in X, then also 

dX , dY ^ . 

^+^=0 (52.) 

(the differentiation with respect to p being in each case performed only so far as p 
appears explicitly in the function). 

The value of Y is given by the equation 

Y=-(X)-|-(«i)yi+(*2)3/2+- + W3/«, • (53.) 

where the brackets indicate that a^i ...a7„ are supposed to be expressed in terms of 

t/, ...3/„ (arts. 2, 3.). 
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Theorem II.— Suppose the function X to contain explicitly, besides the ?^ variables 
Xi...x„, another variable t, and also n constants a^, a^,...a^\ and in addition to the 
equations (50.), let the follovnng be assumed: 

l.ar^^>'-d^=''-^ (54.) 

where b^, ... b^ are n other constants ; so that, by virtue of the 2m equations (50.), (54.), 
the 2m variables x-^...Xn,y.^...y^, may be considered as functions of the In constants 
a,, ... a„, *!,...&„, and t. Then if from the equations (50.), (54.), and their total differ- 
ential coefficients with respect to t, the 2ra constants be eliminated, there will result 
the following In simultaneous differential equations of the first order ; viz. — 

^*=%;'2/i=-^. (55.) 

where Z is a function of .Ti ...x^, y^, ...y„ (which will in general also contain t expli- 
citly), and is given by the equation 

Z=-(f) . (56.) 

In this equation -T7 represents the partial differential coefficient of X taken with 

respect to t so far as t appears explicitly in the original expression for X in terms of 
a;i...a;„, a,...a„ and #; and the brackets indicate that ai, ...a„ are afterwards to be 
expressed in terms of the variables by means of the equations (50.), (arts. 5, 6.) 

Theorem III. — Let the supposition that the 2n variables Xi...x„,y^...y^ are expressed 
in terms of the 2n constants and t, be called Hypothesis I. ; and the converse suppo- 
sition that «!...«„, bi ... bn are expressed in terms of the 2n variables and t, Hypothe- 
sis II. ; then will the following relations subsist : 



dxi dbj dwi duj 

dttj dyi dbj dyi 



dyi dbj dyi da^ 

dttj dxi dbj dx; 



7 



(57.) 



(In each of these equations the first member refers to Hyp. I., and the second to 
Hyp. II. ; and since there is no connexion between the indices of the variables and 
those of the constants, the case of i=j has no peculiarity.) 

Theorem IV. — Let the symbol [/>, q] be an abbreviation for the expression 



« /dp dq dp dq\ 
'\dyi doc~dwi dyi J 



(where p, q are any functions of the 2n variables, which may also contain any other 
quantitieis explicitly ; and the differentiations are performed only so far as a;i,&c.,t/i,&c. 
appear explicitly in p, q); then if «!,...«„, b^, ...b„he expressed {Hyp. II.) in terms of 
the 2» variables and t, the following equations subsist identically : 

[«;, £i] = -[&j, aj = l, [fii, b^=z[ai, a^ = [bi, b;\—0 . . . (58.) 

2 s 2 
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{i being different fromj); and obviously in all cases 

[P> ?] = - [?' ?]> and [p, p]=0 (art. 9.). 
Theorem y. — If u, v be either (1) any two functions whatever of the 2n constants 
fti, &c., bi, &c., or (2) any two functions whatever of the 2« variables a?„ &c., «/„ &c. 
(which may in either case also contain t explicitly)*, then 

y(du dv _du dv^ y j du do du dv~\ CfiQ \ 

\dyidxi dxi dyi] \dai dhi dhi dai] 

(When u, v represent functions of the constants, the differential coefficients in the 
first member of this equation refer to Hyp. IL; and, when functions of the variables, 
those in the second member refer to Hyf. I.) (art. 10.). 

Tlieorem VI. — Let a?„ ... x^, y^, ... y^ be 2w variables concerning- which no supposi- 
tion is made except that they are connected by n equations of the form 

ai—(pi{os„x^,....x„y„y^,...y^) (a.) 

(where the functions on the right are only subject to the condition that the n equa- 
tions (a.) shall be algebraically sufficient to determine 3/1, ...?/„ in terms oi x^,...x„, a^, 
he, and may contain explicitly any other quantities besides x^, &c., y^, &c.). 

Then, if by means of the equations (a.) the n variables 3/,, 3/3, ...«/„ be expressed as 

functions of x„ x^, &c., a„ &c. ; in order that the — s — conditions 

Ayj._dyi 

may subsist identically, it is necessary and sufficient that each of the -^ — - expres- 
sions [tti, aj vanish identically. 

Theorem VII. — Let Z be any function whatever of 2?z variables x^... x„, 3/, ...y„, and 
t. If of the system of 2n simultaneous differential equations of the first order 

. dZ , dZ ,, , 

^^=^,'3/.= -^ (I.) 

there be given n integrals involving n arbitrary constants «!, a^, ... a„, so that each of 
these constants may be expressed as a function pf the variables x^, &c., y'l, &c. (with 

or without t) ; then if the — ^ — conditions [a;, a[\=0 subsist identically, the remaining 

n integrals may be found, as follows. By means of the n given integrals let the n 
variables y^ ...y^ be expressed in terms of i^-i, &c., a^, &c. ; and let (Z) be what Z 
becomes when y, ...y^ are thus expressed. These values oiy^, 3/2 ...3/^ and — (Z), will 
be the partial differential coefficients with respect to Xy, x^, ...x„ and t, of one and the 
same function ; call this function X, then, since its partial differential coefficients are 

* It -was inadvertently stated in art. 10, that u, v must not contain t explicitly. But it is evident that no 
such limitation is implied in the demonstration of the theorem. The preceding theorem is obviously a particular 
case of this ; namely, the case in which u=cti, v—bj. 
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all given (by the equations ^=3/^, -^= — (Z)j, X may be found by simple inte- 
gration, and is therefore to be considered a given function of a?„ ...x„, a„ ... a„ and t. 
The remaining n integrals are then given by the n equations 

bi ... b„ being n new arbitrary constants. 

[On the relation between this theorem and the theories of Sir W, R. Hamilton 
and Jacobi, see arts. 15-20.] 

60. Other results established in the former part will be referred to as occasion may 
require. To the theorems enunciated in the preceding article, the following may 
now be added. 

Returning to the equations (50.), (54.), (55.), we may observe, that if, in thejirsi 
members of (55.), Xi, i/i be supposed expressed in terms of «!, &c., b^, &c. and t, then 

^, -§l may be written instead of *•, ?/• ; since on this hypothesis tlie total differential 

coefficients of Xi, y^ are obtained by differentiating with respect to t as it appears 
explicitly. We have therefore 

doci dTi diji «?Z 

dt dyi dt dxi 

where the first members refer to Hyp. I., and the second to Hyp. II. But since the 
equations (50.), (54.) involve «, h, exactly in the same way as they involve x, y, it is 
obvious that the same reasoning which leads to the equations just written, would 
lead, mutatis mutandis, to the following, which may be considered as an addition to 
the system of equations (57.) (Theorem III.): 

dai dTi dbi dTt /'fiO V 

dt dbi dt dai ' 

In these equations, a^, bi in thej^r*^ members are supposed to be expressed in terms 
of the variables {Hyp. II.), whilst in the second members .^i, &c., y^, &c. are supposed 

to be expressed in terms of the constants and t {Hyp. I.). As before, Z= — -^j but 

in (60.) Z is differently expressed, being what the Z of (55.) becomes when «„ &c., 
3/1, &c. are expressed according to Hyp. I. 

It is to be remembered that all consequences deduced from the form of the system 
(50.), (54.) belong to the system of equations, obtained as in Theorem VII., which 
express the solution of the differential equations (I.). Such a solution will be called, 
as before, a normal solution; and the system of equations obtained by expressing 
fl„ &c., &i, &c. in terms of the variables and t, will be called a system oi normal inte- 
grals. (See art. 20, and the note to art. 29.) 

51. Let «„ &c., 5i, &c. be called, as before, elements. If then c be any function of 
the elements, when the latter are expressed in terms of the variables and t 
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{Hyp. II.), 6" becomes also a function of the same ; and we have 

dc ^ /rfc dai.dc db\ y / dc dZ dc dZ\ /■«! 1 

dt '\dai dt dbi dt / \dai dbi dbi daj 

(see the last article). But, by Theorem V., this becomes 

S=[^'Z]- • • • m 

It is worth observing that both this equation and (60.) might have been obtained 

dc 

indirectly as follows. Since c is constant, we have c'=0 ; that is, -f.-{-[Z, c]=0 (see 

(32.), art. 22.) ; this gives (62.), since [Z, c] = — [c, Z], and again, by Theorem V., is 
changed into (61.); and if, in the latter, we put successively c=aj, c=bj, we obtain 
the system (60.). 

Section V. — On the F^ariation of Elements. 

52. The following general problem includes, I believe, all the cases which occur in 
practice. Let Pi, ... P^, Qi, .... Q„ be any functions whatever of the 2n variables 
Xi, ... x„, y^, ... y^ and t. It is required to express the 2n integrals of the system of 
2n simultaneous differential equations of the first order 

x\=V„y[=Q, (63.) 

in the same form as the integrals (supposed given) of the canonical system 

^^=^/ y^=-d^: (^-^ 

by substituting functions of t for the constant elements of the latter system. 

Suppose a normal solution (see end of art. 50.) of the system (I.) to be employed. 
The elements a^, h^ represent the same functions of x^, &c., y^, &c. and t as before, but 
are now variable ; consequently we have 

, dai fdai , dai ,1 da,- , ^ Irt da, , ^ dat] 

with a similar expression for Z>-. But, by equations (57.) and (60.), these are imme- 
diately transformed into the following : 



■db,-^^^l^-dbr 'db, 



b'=——- 



■' [ ^ dai ■' dttij } 



(E.) 



where Z, Q,-, P^, Xj, y^ in the second members are supposed to be expressed {Hyp. I.) 
in terms of the elements and t. Thus the system (63.) is transformed into a system 
involving the new variables a^, bi, instead of the original variables x^, yi. 

53. If, instead of employing a set of normal integrals of the pattern system (I.), we 
take any complete set of integrals c., c^, ... c^„ then Cj, &c. may be considered as 
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functions of «„ &c., and again, through them, of the variables. We have then 

and if in this equation the values of a';, &c. be introduced from the formula (E.) of 
the last article, the following expression results : 

(in which the symbol {-p, q) is used to denote 



y I dp dq _ dp 
idbu dajt da, 



■k dbk] 



so that by (59.) (Theorem V.) we have {p, q)=i — \_p, q'\ ; but in {p, q] p and q are 
considered as functions of a^, &c., &„ &c., whilst in [/>, q'] they are considered as 
functions of a?„ &c., y„ &c.). Now, considering jo, q as functions of Cj, &c., and through 
these, of «„ &c., we have (by reasoning exactly similar to that employed in deducing 
equation (24.), art. 9.) 

(,,„=.(,.„.,, (1 4-* I)) 

(the summation referring to all binary combinations of the indices a, |8). Hence we 
have, putting 9— c,, 

{;?, cj = 2„({c„, cJ^J, (64.) 

and consequently the above expression for q becomes 

c;.= {Z,cJ+S.2,({c,,cj(Q,|-P,|)), ...... (F.) 

an equation, which is easily seen to become identical with (E.), art. 52, when c^.-.c^,, 
represent ai...a„, &,...&„. 

54. The simplest case is that in which the system of equations (63.), whose inte- 
grals are sought, are of the canonical form ; that is, where 

dyi ' dxi 

W being a given function of the variables (with or without t). In this case the 

formula (E.) becomes 

dZ dW 



' dbi dbi 



' dui dtti - 



(65.) 



whilst (F.) is easily found to be reducible, by the help of (64.), to either of the fol- 
lowing forms : 

c;={Z,cJ-{W, c,} (66.) 
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If we put W=Z+fl, so that Cl may be called the " disturbing function," the above 
formulse become 

<=-f.' *''=S (•=«•> 

c;=S,((c„ cjg) (69.) 

On the first of these forms see the note to art. 38. With respect to the form (69.), 
if we put for {C(, cj its equivalent — [q, cJ, or [c„ c,] (see Theorem V. art. 49.), we 
obtain the well-known expression 



, „ /' -,dD,\ 



The difference between this last form and (69.) consists in this ; that in the latter the 
coefficients [c„, cj are obtained from the expressions for c„ c^, &c. in terms of the 
variables; whereas in (69.) the coefficients {c^, c„} are similarly obtained from the 
expressions for c„ &c. in terms of the normal elements a^, &c., b^, &c.* ; and when a 
normal solution of the undisturbed problem has been obtained, the latter process will 
generally be found much more convenient than the former, since the elements c,, &c. 
will usually be much simpler functions of the normal elements than of the variables. 

55. In illustration of this, it will be worth while to deduce the expressions for the 
variations of the ordinary elliptic elements of a planet's orbit from those of the normal 
elements given in art. 30, 

Let a and e be the semiaxis major and excentricity, / the inclination of the orbit to 
a fixed ecliptic, v the longitude of the node, w the longitude of the perihelion, wi+(g) 
the mean longitude of the planet ; longitudes being reckoned in the plane of the 
ecliptic (from a fixed origin) as far as the node, and then on the plane of the orbit . As 

usual, « stands for —• Also let n^+(2)=l «rf^+s, so that g'=(g)'+j?n'. 
If, then, we call the six normal elements Kj, a^, Wj, /3i, jS^, /Bg, we have (see art, 30.) 

Kg = m^^a{ 1 — e^) . cos /, ^^=y; 

from which, conversely, 

(0=^.+ft-?#A. 



'*— 2ai 


5 


1 «2 — 


2«i«l 


1 — e^=: 


■mV' 


cos< = 


«3 

■ — 5 



•JCj, Oal — iijl -— — _. -— . 1. 

\abf daj daj dojj 
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From these expressions the values of {a, e}, {a, /}, &c. are found with the greatest 
simplicity, and the results are 



m(j,{a, (£)}=2wa^ m(j,{{B), e}=^?^-^^^^(l-x/l-e^). 



TO;M,{w, e} =-—^ -, m^{{t), t} =z-~^= tan -, 

muj{is, t) ■=. ,- — -i.tan-i mujh,t}=:-. . , ? 

the rest all vanishing. Hence, observing that if R be taken in its usual signification 
we have fi =— R, we obtain* 

K0=-2«a^^+—^—(l-yi-e^) ^+^7137^ tan--, 
, -«« f 1 dR , . I (dR , dR\\ 



[/jV 



na dR 

sintA/l — e^ «?» 



3 



in which we may, as usual, put s for (s), provided that in forming the term -^, nt be 

exempt from differentiation with respect to a. 

56. A comparison of the above JDrocess with that by which the corresponding 

* If we consider R as a function oip, q instead of i, v, where |)=tan ( cos t>, gi=tan i sin v, we find 

dB, / dR . dR\ 

^=sec%(^cosv^+sinv^j 



cosy-^ sin v 

dv ' 

and consequently 



dp dq , 

= tan I ( c 



dn . / dR . dR\ 

•— -= tan I ( cos y — sin v-r~], 

dv V dq dp J 



, — M«!(seci)^r dR,. I (dR ,dR\\ 

/»»'= — , . <^sec( F-tan-cosvf— — + -7- \ 

'^ t/l-.e'i\ dq 2 \d{s) dwj j 

,_na(seci)^( dR ( . /dRdR\\ 



The formulae will then agree with those of the M^canique Celeste (Supplement to vol. iii. p. 360, ed. 1844), if 
we allow for the diiferent mode of measuring longitudes, and neglect, as Laplace does, terms of the second 

order -Urith respect to i andi^. (Lapiace uses R with the opposite sign.) Those in the text agree (allowing 

for notation) with the expressions given by Professor Hansen, Astr. Nachr. No. 166, art. 3, equations (2). 
MDCCCLV. 2 T 
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expressions are obtained by Pontecoulant*, will show the convenience of using the 
coefficients {Cj, c,} instead of [c^, c,] (in Ponte'coulant's notation (c^, c,)). 

[It will be observed that the formulae for (s)', is', t' at the end of the last article, 
do not agree with those of Pontecoulant (p. 330) for the variations of the corre- 
sponding quantities a, &>, <p. The reason of this is as follows : — In Pontecoulant's 
notation <p expresses the same as / in this paper, and a the same as v. But co (the 
longitude of the perihelion) is not the same as w; the former being measured entirely 
in the plane of the orbit from a radius vector, ^xed in that jdane-f, and assumed as 
the origin of longitudes. Consequently s, in Pontecoulant (which we will call g, for 
distinction), is not the same as (e) in the present paper. In fact, if we equate the 
expressions for the mean anomaly in the two notations, we have 

also it is evident that if we put /3 for the angle between the node and the origin from 
which 6) is measured, we have (1(5= — cos tdv, and zir=i'-j-/3-|-iii', so that 

d'^=d&)-\-{l — cos i)dv. 

If' then it were allowable to consider R as capable of being expressed as a function of 
0) and Si instead of •m and (s), and if we represented by (R) the expression for R so 
transformed, we should have 

d^dr.+^^d(e) + kc.=^-^da>+-^ds,+kc.; 

and if, in the two first terms, we put for dwr and d{s) the values d'^sy = dai-}- {I— cos i)dii, 

d{s)=dSi-\-(l— cos t)dv, and compare the two expressions, we find 

^_rf(R) dR_d{R) 

d'U! dm ' d{e) ds^ - , 

dR , . ./dR , dR\ 

„+(l_cosO(^^+^j: 

These relations, together with the equation 

J</ = ii)'+(l— COS/)/, 

are easily seen to render the expressions at the end of art. 55 identical with those of 
Pontecoulant ; in fact, it is by an equivalent transformation that the latter are 
finally obtained by that author from the correct expressions in p. 328. But it is to 
be observed that this proceeding is founded upon a false assumption ; for it is im- 
possible to express R as a function of a, e, t, v, s„ u, as is obvious from the considera- 
tion that R, in its original form, is not a function of (g)— w merely, but also of (s) ; 
whilst (s) is not expressible as a function of the new elements, as is shown by the equa- 
tion </(s)=rfe^+(l—cos/)6?j':}:. It would be out of place to enter further into this sub- 

* Th^orie Anal, du Systfeme du Monde, tome i. pp. 316-330. 

f On the meaning of this expression, see below, art. 73. 

X It -would be a work of some trouble to trace accurately the process by which Laplace arrives at the for- 



d{R) 
dv 



V 
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ject here, especially as some of the most important principles involved in it have 
been discussed elsewhere *. See also Appendix B.] 

57. Returning to the expression (69.), art. 54, it may be observed that the coeffi- 
cients {Ci, Cj} are to be expressed in terms of Cj, c^, &c., and this involves no difficulty 
when each of the two sets of elements Cj, &c., a^, &c. can be expressed in terms of the 
other explicitly, as was the case in the example just discussed. Suppose, however, 
that the normal set a^, &c., &j, &c. are given in terms of the set Cj, &c., but that it is 
impracticable or inconvenient to obtain the converse equations expressing the latter 
in terms of the former. In this case we may proceed as follows. 

Adopting the notation of art. l-f", and putting/, g for any two of the set Cj, c^, &c., 
we have 

suppose this equation written at length, and then, after multiplying by ' ■','^' , let 

each side be summed with respect to all binary combinations/,^. The result is (see 
art. 1, equation (4.)), 

('j|ff-{/.S}) = > (70.) 

(the summation referring to the combinations / g). Again, if the former equation 
be nmltiplied by ^,J ! , where p, q represent any two of the normal elements, a„ &c. 
fcj, &c., except a conjugate pair, and the sum be taken as before, we have 

The two formulae (70.), (71.) give n{2n—l) linear equations for determining the 
n{2n—l) unknown quantities {f,g}; the coefficients of the latter being all given 
functions of c^, &c. But such cases will hardly occur in practice. (With respect to 
the form of the above system of linear equations, it is easy to show that the complete 
determinant of the coefficients is =1.) 

58. The integration of the formulae (65.), art. 53, would give the means of ex- 
pressing the solution of the system 

,_dW ,___dW 
' dyi^ y* dxi 

mulse alluded to in a preceding note, as the various steps of it are to be found in different places, the notation 

is somewhat inconsistent, and the results do not profess to be rigorous. My impression is, however, that 

Laplace nowhere commits the fallacy of assuming (for example) that R is a function of r, v, z, or r, v, s (see 

vol. i. p. 295), where v is the angle described by the radius vector on the varying plane of the orbit. 

* See Jacobi's two letters to Professor Hansen in Cebllb's Journal, vol. xlii. 

, . . d(u, v) , , . ^. J. du dv du dv 

t t. e. using -) (■ as an abbreviation for — - — — • 

d{x, y) dw ay dy ax 

2 t2 
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in the Jbrm of a normal solution of any other similar system 

'—^ ,__^ 

' dy' "' dx' 

which may be chosen as the pattern. 

In the most usual examples the function to be chosen for Z is naturally suggested 
by the circumstance, that W presents itself under the form of the sum of two func- 
tions Z+n, of which the former, taken alone, gives an integrable system. But this 
is not necessarily the case; and it is worth while to observe that the formulte (65.) 
take a simple and remarkable form whatever Z may be, provided that it be a func- 
tion not containing t explicitly. For then, assuming the " integral of vis viva," Z=.h, 
as one of the normal integrals of the pattern system*, the element conjugate to h is 
7 (the constant added to t) ; and observing that Z, in (65.), being expressed in terms 

dZ 
of the elements, reduces itself simply to h, we shall have ^=1, whilst the differential 

coefficients of Z with respect to all the other elements vanish; so that, if we put 
«i5 •••««-» bi,---K-i for the remaining elements, the system (65.) takes the following 

dW y_dW j 



(72.) 



"'■■ dbi' "i— dai j 

This, in dynamics, gives the process to be used in the following problem : " To express 
the solution of any dynamical problem in the form of the solution of any other [involving 
the same number of variables) in which the principle of vis viva subsists^ 

59. As an example of the above process we may apply it to determine the motion 
of a simple free pendulum (not taking into account the earth's rotation). 

Let / be the length of the pendulum, and let the mass of the material point m 
placed at its extremity be represented by unity. Also let x, y, z be the rectangular 
coordinates of m, the origin being at the position of rest of m, and the axis of z 
directed vertically upwards. The equation to the sphere described by m is 

and the force-function U is —gz. 

Hence if we take, as the two independent coordinates, the radius vector § of the 
projection of m on the plane of xy, and the angle 6 between § and the axis of x, we 
shall have for the differential equations of motion, 

(A.) 



du' ' dv 



, dW , dW 
See art. 19 (where hi in equation (29.) is a misprint for Jj). 
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where u, v, are the variables conjugate respectively to f, S, and defined by the equations 

and W is T— U expressed in terms off, d, u, v. 
Now x=§ cos d, t/=§ sin 6, z=l—^/f—f ; hence 

from which the following expression for W is easily obtained : 

-w=l(^^u^+j^+g{i-Vi^). ...... (W.) 

Now let us take as a model for the solution of the above system, a set of normal 
integrals (in polar coordinates) of the system 

x"+w^x=0, y"+n''i/=0, (B.) 

where n^ =7- In this system we have U=— o'^Y; and proceeding exactly as in 
art. 27, we obtain the following results : the two integrals of vis viva and of areas are 

h=l[u^+f.+nY)] ^._^ 

c=v J 

these are to be solved for u, vi and then V is to be obtained from the equation 
cfV =udg+vdd. This gives 



V=zc6+(dpl2h-nY-i]* ; 



r 

and the remaining integrals are given by the equations 

dY , , dY 

r and w being the elements conjugate respectively to h and c. Performing the dif- 
ferentiations ^^r*^, and taking the integrals in the second term so as to yanish with 

the expression {^h-n^-^^ (seeAppendi. A.), we find easily the final equations 

«V=^+\//i'— wV.cos 2w(^+t)1 

cV=A-y^^-/iV.cos 2(^-^)1' 

in which w is the angle between the axis of x and a (distant) apse, and — r is the 
time of passage through that apse. The four equations (i.), (ii.) comprise a com- 
plete normal solution of the equations (B.). The last is the polar equation to the 
elliptic orbit ; and if we call a, h the semiaxes of the ellipse, we have 

c=nab, h=nr — - — 
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60, The solution of the system (A.) of the last article will now be expressed by the 
same equations (1.), (ii.), if the elements A, c, r, m be variables defined by the system 
(see art. 58.) 

ar ah 

d'sr do 

where W is to be obtained by substituting in the expression (W.), art. 59, the 
values off, 0, n, v in terms of the elements and t, derived from equations (i.), (ii.). 
The result of this substitution is 

+^3? cos 4«(^+r) H-n^/(/- x/F=7), 

in which the value of f^ in the last term must be understood to be substituted from 
the first of equations (ii.). if we call p the angle between the pendulum and the 
vertical, we shall have evidently 

nH{l-^/K^) - nH\ 1 - cos ?>), 



and the differential coefficient of this term with respect to any constant h involved 

in the value of p will be -^ A^. Obi 

" 2 cos if) dk 

sions for the variations of the elements : 



in the value off will be -^ fjj- Observing this, we obtain the following expres- 

^ COS ^ diC 



m- 



A'=-«(secf-l)yA^— wV.sin2n(if+r) + ^— |i^sin4n(;+r) 



nC' 



^"-2lr^+^(«^^^-^)(^+:^l^-^««2«(^+^))+-^,cos4«(^+.) 

c' = 

-'=1 -^(sec 9-1)7^, cos 2n{nr) -^, cos An{t+r). 

The third of these equations gives ab=constant ; hence, by means of the equations 
at the end of art. 59, the following expressions are easily deduced : 

-= —j—^\— (sec p- 1) sin 2n{t+r)-\—jr~ sin An{t-{-T) \ 

r'=-^p-(- 1 + cos 4n(#+T))+2(sec ?!- l)n +^5Z:p cos 2??.(/+t) j 

, , fl— cos4«(^ + r) , , . cos2w(^ + t)"1 

n/—nab]^ 2/« - — (sec^— 1) ^^^^^ ' V- 

These equations are rigorous, and in general not easier to integrate than the original 
system of which they are a transformation ; but they may be integrated approxi- 
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rnately on particular hypotheses. For instance, if the pendulum never deviates much 
from the vertical, % is always small, and sec (p— 1= -|™ nearly ; introducing this value, 

and substituting for f, we have for the non-periodic parts of the above expressions 

b' 



a 


-b=^ 




r'=- 


a' + b^ 


/ 3 nab 



Hence the axes of the mean ellipse are constant ; also if T be the time of describing 

2t 1 

the ellipse, we shall have approximately ^=—'^3775 whence it follows that the mo- 
tion of the apse during this period will be -^- -tt— t\' ^"^ "^"^"^ nearly. 

This agrees with the statement of the Astronomer Royal*. The above approxima- 
tions would cease to be sufficiently accurate if b were very nearly equal to a, or the 
motion very nearly circular ; but they apply to any other case, and in particular to 
that in which b is very small, or the motion nearly rectilinear. In the case of nearly 
circular motion, it would be necessary to develope sec p in a series of powers of 
a^—b^^ and the results would be applicable whether p were small or not. But I shall 
not pursue this subject further here. 

Section VI. — On the Transformation of F^ariables. 

61. The method of the variation of elements, theoretically considered, consists 
merely in a transformation of variables of a particular kind ; that kind namely, 
which leads to a new system of differential equations belonging to the same general 
class as the original system. But practically, the choice of variables is determined 
by the well-known considerations from which the method derives its name. 

It is the object of the present section to consider the general class of trans- 
formations of which the method in question is a particular, and not the only useful 
case. 

62. Definition of Normal Transformations. 

Let li, I2, ... Lj 'Jd ^'2, ••• ^n be new variables connected with the original variables 
Xi, &C.3/1, &c. by 2n equations (which may also involve t explicitly), such that each 
variable of either set may be considered as a function of the variables of the other 
set (with or without t). Let P be any function of |j, I2, ... |„ y^, y^, ...y^ and t ; then 
if the equations connecting the old and new variables can be put in the form 

I propose to call the transformation normal. 

* Proceedings of the Royal Astronomical Society, vol. xi. p. 160. 
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[It follows from Theorem I. art. 49, that the system (73.) is equivalent to each of 

the following': 

dQ_ dQ,_ 

in which Q= — P-\-'2i{Xiyi), and is expressed in terms of x„ ... x„, |i, ... i„; or 

dR_^ dR_ 

in which R= — P+Sj(liJj,), and is expressed in terms ofy„ ...«/„, ^ji, ... jj„; or lastly, 

dS__ ■ d^__ 

dru"^'' dx~y^^ 

in which S= — P+2i(a7iyi4-S.%)5 and is expressed in terms of a',, ... sc^, rii, ... ^„. 

Any one of these forms might be used ; but I shall employ the form (73.) for 
reasons of convenience.] 

63. Inasmuch as the equations (73.) of the last article are of the same general 
form as the system (54.), art. 49, all the conclusions deduced from that form will 
subsist, mutatis mutandis ; so that we may apply the Theorems (II.), (HI.), (IV.), (V.), 
art. 49, by merely changing X into P, and 

Xi, ... x„,i/„....i/„, a„ ... a„, b„ ... b„ respectively into 



"n » 



observing that instead of 4> y't we must now write -^, J?*. We thus obtain the fol- 
lowing relations : 

dt dm' dt d^i' • ^'"^-^ 

where "*F is a function of fi, &c., ;;„ &c. and t, defined by the equation 



•qr 



= -(f). (7^.) 

the brackets indicating that the expressions for y^, ... t/„ in terms of the new variables 

dP 
li, &c., ;?!, &c., are to be substituted in -^ after the differentiation with respect to t; 

which is performed so far as t appears explicitly in the original expression for P as a 
function of |, ... |„ «/i ... y„, and t. (See Theorem II.) 
We have also the system 

d^i , dXf d^i d-i/j' 

dyj drii dxj dy\i 



drii dwj- drii di/j 



dyj d^i dxj d^i- 



(76.) 



: * For in the original theorems Xf is the same thing as the differential coefficient of ,Vi taken with respect to t, 
as t appears explicitly in the expression for iCi in terms of «j, &c., bp &c. and t ; the analogous quantity in the 
present case is therefore the diflferential coefficient of f « taken with respect to t, as t appears explicitly in the 
expression for ^i in terms of x^, &.C., y^, &c. and t. But this must not now be denoted by ^^, inasmuch as 
iTi, &c., yj, &c. are themselves afterwards to be considered as functions of t. 
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(see Theorem III.). And if p, q be any two functions of the variables (with or with- 
out f), then 

'\dyidXi dxidyj \drti d^i d^idrnj' ' ' ' ' ' ^' ' '' 

where p and q in the first member are supposed to be expressed in terms of vTi, &c., 
ifi, &C.J and, in the second, in terms of |„ &c., ??„ &c. In other words, the value of 
[p, q\ is the same, whether it be obtained from the expressions for p, q in terms of 
the original variables, or by an analogous process from their expressions in terms 
of the new. 
Particular cases of (770 are the relations 

[i,,s,,] = -l, [l„fj = fe»;j = [l-,|]=0 (78.) 

(See Theorems IV., V.) 

64. The relations (74.), {7Q.), {77-), (78.) of the last article depend solely upon the 
form of the equations (73.), art. 62, which connect the new variables with the old ; 
and are independent of any supposition as to the equations which may determine 
either set of variables as functions of t. Let us now, however, introduce the suppo- 
sition that the original variables Xi,.,. x„, t/,, ...y„ are determined as functions of t by 
the system of differential equations, 

' dZ I dZ /T \ 

^'=d^: ^^=-^, ('•) 

The relations just established enable us immediately to transform this system into 
another involving the new variables instead of the old ; for we have 



s=f +s<|%+|3<) ■■ 



now §=:^ (see (74.), art. 63); 

and if in the remaining term we substitute for x'j, y\ their values from (I.), and for 

$i, $i their values ^, — ^, it becomes 
dxj dyj dr\i aijj 



^ /dZ dxj.dZ dyj\ 
\dxj di\i dyj drtJ' 



dZ 



which is equivalent to -j-} if Z be supposed expressed in terms of the new variables. 

We have then 

and, exactly in the same way, 



We have then l<=-r-+T-' 

dni dr)i 



''* d^i d^i 
This result may be stated in the form of the following Theorem VIII*. If the system 

* This theorem, in its general form, is, to the best of my knowledge, new. But that case of it in which P 
does not contain t explicitly has already been proved in a different way by M. Desbovbs, who has, by means 
MDCCCLV. 2 U 
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of differential equations (I.) be transformed by the introduction of new variables 
§1, ... |„, ;?!, ...;?„, connected with the original variables x^, ... x^, y^, ...yn, by the equa- 
tions ■j-=Si, ■yy=^i, where P is any function of li, ... |„, ?/„ ..-yn, which may also con- 
tain t explicitly, then the transformed equations are 

in which O is defined by the equation 

at 
and is to be expressed in terms of the new variables. (The substitution of the new 
variables in -^ is to be made after the differentiation. See art. 63.). 

Corollary. — If P do not contain t explicitly, ^=0 and ^=Z ; so that in this case 

the transformation is effected merely by expressing Z in terms of the new variables. 

65. It follows from (77 .)j art 63, that \if,g be any two integrals of the system (I.), 
the value of [/, g] is the same whether it be derived from these integrals in their 
original form, or similarly obtained from the same integrals after transformation by 
the introduction of the new variables. And consequently if n integrals a„ aa, ... a„ o/" 

the original system be given, which satisfy the-^ — - conditions [a^, aj]=0, they will 

continue, after a normal transformation, to satisfy the analogous conditions, so that the 
method of finding the remaining integrals given in Theorem VII. art. 49, will also 
continue to be applicable. We had an instance of this in the case of the problem of 
central forces (art. 27.) j where the above conditions were found to subsist after the 
transformation from rectangular to polar coordinates. (It will be shown presently 
that every transformation of coordinates is a normal transformation.) 

66. It was shown in Part I. (art. 18.), that if W be any function of x„x^, ... x^, 
x[, 4, ...x'„ (which may also contain t explicitly), the system of n differential equa- 
tions of the second order 

\dx'i / doci \ ') 

may be changed into a system of 2n equations of the first order of the form (I.), 

of it, deduced Jacobi's form of the method of the Variation of Elements (namely, the equations (68.), art. 54), 
from the similar form of Lagrange, in which the elements are the initial values of the variables. It will appear 
in. the sequel that the extension to the case in which P may contain t, is of importance. If the expression were 
not already appropriated, I should have proposed definitively to call P the " modulus of transformation j" and I 
shall use this term provisionally in the present paper, not being able to suggest a tolerable substitute. After all, 
as the word " modulus " itself is used without confusion in very different senses according to the subject matter, 
there is, perhaps, no reason why a similar liberty should not be allowed in the use of the proposed expression. 
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art. 64, by putting 3/j=-^5 and taking Z a function of x^, &c., 3/1, &c. and t, defined 

i 

by the equation Z=— W+Sj(«-yi), 

in which x\, ...x'„ are to be expressed in terms of 3/1, &c., *i, &c.*. 
Conversely, a system of the form (I.) being given, it may be changed into a system 

of the form (80.) as follows: by means of the equations a;.=— , let y,, ...?/„ be ex- 

pressed in terms ofx[, &c., Xi, &c. ; it follows from Theorem I. (art. 49.) that we shall 

have i/,=^ . • . (a.) 

and ':r—~'7~' (b.) 

where W is a function of x„ &c., x\, &c. defined by equation 

W=-Z+%{x[y^, 
in which 3/1, ...y^ are to be expressed in terms of x\, &c., x„ &c. The n equations 
y\——^ are then changed by (a.) and (b.) into the form (80.). 

On the Transformation of Coordinates. 

67. It has been seen in the preceding article, that we can always change the system 
of 2n equations of the first order of the form (I.), art. 64, into a system of n equa- 
tions of the second order of the form (80.). In this latter form the equations of 
dynamics naturally present themselves. 

Now in the ease of the dynamical equations, Xi,x^,...x^ are the independent 
coordinates of the system (the word coordinates being taken in its most general sense), 
and when the equations are to be changed into the form (I.), the additional variables 

3/1, ...y„ are defined by the equations -^=3/^. In this case a transformation of coordi- 

nates, in the most general sense, consists in taking n new variables li, I2, ... |„, con- 
nected with the original coordinates Xi, ...x„ by n equations, which may also involve / 
explicitly. It is a well-known theorem, that the transformation of the equations (80.) 
is effected merely by expressing W in terms of the new coordinates |„ ...S„, and their 
differential coefficients I'l, ... C instead of the old ; so that the new equations are 



/dW\< dW 

{<)=^-- <«'•> 



the proof of this theorem does not depend upon the form of the function W ; and we 
know also (see arts. 18 & 66.), that whatever be the form of W, these new equa- 
tions may be again transformed to a system of the form (I.), by taking n additional 

* This theorem is a generalization of Sir W. R. Hamilton's transformation of the Dynamical Equations, 
See Part I. art. 18. 

2 u 2 
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variables fii, ...}]„, defined by the equations 

dW 

[It is to be observed that this last transformation is not, in general, equivalent 
to expressing the original Z in terras of the new variables ; for the original Z 
is — W+2(x^i) (art. 66.), and the analogous expression derived from (81.) is 
— W+2(|j;?i), which is not, in general, equivalent to the former. It will be seen 
presently that the two expressions are equivalent when the equations connecting 
x„ ... x„ with ?i, ... |„ do not involve t explicitly.] 

68. Following the analogy of the dynamical equations, I shall adopt the following 
as the 

Definition of a transformation of coordinates. 

The original equations (I.), art. 64, having been changed into the form (80.), art. 66, 
let li, laj ••• t» be n new variables connected with the n variables x^, x^,... *„ by n equa- 
tions, which may also involve t explicitly ; and let ??„ n2> ...?;„ be h other new variables 

dW 
defined by the equations -^=;?,. (where W has been expressed in terms of §„ &c., S'„ &c.). 

By means of the 2n assumed relations, the 2n original variables x^, ...x^, y^,...y^ can 
be expressed as functions of the 2w new variables li, ...5„, ??i, ...??„. Let this substi- 
tution be called a " transformation of coordinates." 

It has been seen in the preceding article that the original equations (I.) are changed 
by a transformation of coordinates into a system of the same form, which however 
cannot in general be obtained by merely expressing Z in terms of the new variables. 
But we are not at liberty to assume (and it is not generally true) that a change of 
the system (I.) into another of the same form is a normal transformation (art. 62.). It 
has already been stated, however (end of art. 65.), that this is true in the present case ; 
a proposition which I proceed to establish. 

69. Every transformation of coordinates is a normal transformation. 

To prove this theorem, we have to show that every transformation of the kind 
described in the last article is also of the kind defined in art, 62; in other words, 
that it is possible to assign a function P, of f„ ... ^, y„ ...y„ (with or without t), such 
that the given relations between x„ &c., |„ &c., which define the transformation of 
coordinates, shall be equivalent to the system of equations 

dF_ dP_ 

Take P=(^i)3/,+ (^.)3/2+...+(«„)3/„ (82.) 

(the brackets indicating that x^,...x^ are to be expressed in terms of f„ ... |» ; so that 
P is a function of ?„ •••!« «/» —yn, with or without t according as the equations con- 
necting x„...*„ with |„ ...|„ do or do not contain t explicitly). Then P is the func- 
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dV 
tion required. For we have, at once, j-=Xi ; also 

d^i-^' d^i ^y^ d^i +-+y» rff. ' 

now the definition of ra is tii=-^, where 

w=-(z)+%,+%2+...+a?;t^„ 

(expressed in termsof |i, ...|„, |'i, ...|^ (art. 66.)) ; hence, putting Z (without brackets) 
for the original form of Z in terms of «,, &c.,3/i, &c,, and observing that (Z) becomes 
a function of I',, &c., only through i/^, &c., we have 

dW dZd^ dZdy^ _dZd^ 

d^'.—" di/i d^\ ~ %2 d^'~'-- dy„ d^'. 

The two first lines of this expression vanish by virtue of the equations ^i=j-, ; and 
since ^'.'fe) , ^r j.'^^^'^^'-l j_^)?' 

, dx, dix.) 

we have — f=:~7^; 

so that we have, finally, 

_^W_ ^(«i) , d{x^ d{x,^ 

"i-d^'. -y^ dii "1"^^ d'c.i +-+:y« 4, ' 

which is evidently equivalent (see the expression (82.)) to 

the proposition in question is thus established, and may be enunciated as fol- 
lows : — 

70. Theorem IX. — Every transformation of coordinates is a normal transformation, 
of which the modulus* P is a function of^^, ...|„, 7/1, ...«/„ {with or without t) given hy 
the equation P=(^03/:+(^.)3/.+ - + Wy» 

{the brackets indicating that Ki,...x„ are to be expressed in terms of^j,...^„). (See 
arts. 62, 68.) 

This thtiorem will be made more intelligible by applying it to a very simple 
example. 

Let it be proposed, then, to transform from rectangular to polar coordinates the 
differential equations of any dynamical problem referring to the motion of a single 
material point whose mass is m. Let x, y, z be the rectangular coordinates of m, 

* See note on Theorem VIII. art. 64. 
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and u, V, w the variables conjugate to them ; so that, putting T=^m(x'^+y'^+z'^), we 
have u=-T-;,&c.; whence T=2^(u*+v^+w^)j f^nd the equations of motion are x'=-^> 

(17 1 

u'= —7-5 &C.J where Z=2~'("^+v^+'*^^) ~"U, and U is a given function of x, y, z, with 

or without t. 

Now let r, 6, <p be polar coordinates of P, so that 

x=:rsin^cos?), y=r sin ^ simp, z=rcos^. 

Let u, V, w be the variables conjugate to r, $, (p. Then the ordinary process of trans- 
formation would be as follows : — 

(1) to express x', y', z' in terms of r, 6, (p, r', &, (p', and thus transform T into a 
function of the latter quantities ; 

(2) to define u, v, w by the equations 

dT dT dT 

and by means of these relations to express r', 6', (p' in terms of u, v, w, r, 6, <p, so 
that x', y', z', and therefore, finally, T and Z, might be expressed as functions of 
the six new variables. 

Instead of this, let us adopt the method indicated by the theorem at the beginning 
of this article. 

We have then, for the modulus of transformation, 

P=(x)u+(y)v+(z)w, 

in which x, y, z are to be expressed in terms of r, 0,(p; so that the proper form of P is 

P=ur sin 6 cos ip+vr sin 6 sin (p-\-wr cos d, 

dF 
and the equations (corresponding to ^i—^: i^^^- 69.)) which define the new variables 

u,v,w,ave ^p ^p ^p 

^=^' ^=-J-' «^="^- 

These give M=sin ^(ucos?)+vsin^)-|-wcos^, 

v=r cos ^(u cos (p-\- V sin <p) —rw sin 6, 

w= —r sin ^(u simp— v cos p), 

from which the values of u, v, w are easily obtained in terms of the six new variables. 
But in order to eflTect the transformation of T, we have only to square each side of 
these equations, and add them, after dividing the second by r^ and the third by 
(rsin^)^; we thus obtain 
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and the transformed equations of motion are therefore 

'mr'=.u, m^=vr~^, m<p' =w{r mi Oy^ , 



mv'-={r?h.\\S)~^w^co%6-\-m-M-i w'— 



d\J , dU 
d<p' 



The preceding process is not, of course, given for the sake of the result, which may 
easily be verified directly, but in order to illustrate the meaning of the theorem on 
which it depends. It is hardly necessary to add, that if the problem involved the 
independent rectangular coordinates of any number of material points, the trans- 
formation to polar coordinates would be effected in the same way, by merely adding 
to P analogous terms for each point. 

71. Before proceeding further there is an important remark to be made. 

It has been hitherto assumed that the modulus of transformation P was a function 
of no other quantities than the 2n variables f„ ... |„, 3/1, ...i/„ and t. But if every step 
of the demonstrations of the theorems of transformation which have been given (art. 
62, &c.) be examined, it will be seen that they continue to hold good in the following 
more general form. 

Take P=/(ii, I25 • • • L, y^^V'^,— Vn, p, ?, ^, • • •, 0> 

where p, q, r, ... are any functions of any or all the variables, old and new, with or 
without t. 

Let the equations connecting the old and new variables be, as before, 

^=^'' C='''' ^^^-^ 

with the condition that p, q, r, ... are exempt from differentiation informing these equa- 
tions. 
Then take "^— — (r^if)^ ^^^^ the following signification; (1) -^ denotes the 

differential coefficient of P with respect to t, so far as t is contained explicitly, and 

also through the variables in p, q, r, ... ; that is to say, 

d^_d^ dV , ^ , 
dt —dt'^dp'P^dq'9'^'- 

(where p'=-j.-[-^ Xi-\-&c. &c., but this substitution is not to be made at this stage) ; 

(2) (-\if) denotes the result of substituting in the above expression the values of 

y„...y„ in terms of li,...L, r,y,...n„,p,q,r,..Ai-om (83.), so far as —^ contains 3/,,... t/„ 
explicitly {i. e. not involved in p, q, r, ...). 

Lastly, take 0=(Z)+-*-, 

where (Z) denotes the result of substituting in Z the values of the old variables as 
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given in terms of Sj, ... ;?„ ... jo, q, r, ... from (82.). We shall then have 



d^ ' d^ (Q. s. 



vphere is in general a function of 

f„ ... |„, ??„ ... }j„, p, q, r, ... p', q', r', ... and t; 

and the differentiations with respect to |„ ;?; are performed only so far as those vari- 
ables appear explicitly in O. But after these differentiations, we may introduce the 
actual values ofp, ... jo', ... in terms of the variables and their differential coefficients. 
It is obvious that the original variables will not, in general, have been eliminated 
from the system (84.) ; but of course the elimination may be afterwards completed*. 
Similar considerations apply to that particular case of transformation which we have 
called a transformation of coordinates (art. 68.). We have then 

and the relations connecting x^ ,.. x„ with |„ ... |„, may contain p, p, r, ... ; so that 
(^1), &c. are functions of j», q, r, ... as well as of li, ... |„. 

We might have deduced the preceding conclusions from the following simple con- 
sideration. Since p, q, r, ... are actually functions of t, though unknown functions, 
we may imagine them to be known, and to be expressed explicitly in terms of t ; and 
then the case resolves itself into that of art. 62, &c,, so far as the demonstration is 
concerned. But as a doubt might possibly have arisen whether any fallacy was 
involved in the circumstance that p,q,r,... involve (when supposed to be expressed 
in terms of t) the arbitrary constants of the problem itself, it seemed best to refer to 
the original reasoning ; the most important part of which is that contained {mutatis 
mutandis) in art. 6. (For the " mutanda" see the beginning of art. 63.) It is then 
apparent that this circumstance is perfectly immaterial with reference to the con- 
clusions in question, though it may be important in other points of view. 

72. This being premised, we will proceed to an example of transformation more 
interesting than the former, namely, the 

Transformation from fixed to moving axes of coordinates. 

Let x, y, z, u, v, w have the same signification as in art. 70, and let x,y, z, u, v, w 
be the new variables, where x, y, z are rectangular coordinates, referring to a system 
of moving axes of which the origin always coincides with that of the original fixed 
axes of X, y, z. 

Let the direction cosines of the new (moving) axes with respect to the old be 
\, (^0, "0 ; K, ^i, h ; K, /*2, h, thusf- : 

* The final equations in this case will not In general have the canonical form. 

f I do not know who first used this convenient way of indicating the nine direction cosines by a diagram, 
but I first saw it in one of M. Lame's works. 
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y^. 
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^0 


i^x 


f*2 


z 


"o 


h 


"2 



where X„, &c. are functions of t, which may be either given explicitly, or implicitly 
through the variables (see the last article). 

The modulus of transformation P is found (art. 69.) by substituting for the vari- 
ables X, y, z in the expression xu+yv+zw, their values in terms of x,y, z; we have 
therefore 

and then the three equations 

dV dF dF . 

^=^*, -^=v, -^=w, give 



V ■=\^n-\-^{v-\-v{w j-, whence < v =iJij^u-\-[/jiV-\-(^iW 

Lw= ^0^+ "1^+ "iW^ 



(85.) 



Also we have (see Theorem VIII. art. 64.), since P is to be considered to contain t 
explicitly through X„, &c., only, 



dF 

-^ = {>:^x-\-l!,y-\-X^z) u + (|M,>+^;y +^-22;) V +{v<^+hy-\-h^) h^, 

in which expression the values of u, v, w in terms of the new variables are to be sub- 
stituted from (85.). Now if we put ea^, »i, w^ for the angular velocities of the moving 
system of axes about the axes of x, y, z, respectively, so that 

it will be immediately seen that the usual relations between the nine direction cosines 
enable us to put the result of the substitution in the following form : 

\-^j ^=(»)(,{yw—zv)-\-ci)^{zu—xw) -\-6)^{xv—yu). 

The original differential equations 

X'=-^-j u'=— -^5&C 
«u ax 



are then transformed into 
where 

MDCCCLV. 



x=-j-} m'=— -—, &c. (art. 64.), 
dF\ 



*=(Z)-(w) 



2 X 
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Introducing the above value of (-jr), and omitting- the brackets, we obtain for the 
system of transformed equations, 



X 



,' — 



dZ , , dZ , 



, dZ . , dZ . 

, dZ . , dZ , 

aw "^ dz 



(86.) 



in which Z is supposed to be expressed in terras of the new variables. 

73. On the principles of the integration of this, and of transformed systems in 
general, I shall make some remarks hereafter. For the present, the following may 
be observed. If, in the transformation of the last article, we suppose the motion of 
the new axes given, then \, &c., and therefore also w^, a^, u^, are given explicit func- 
tions of t. But if the motion of the new axes is only given by connecting it with the 
motion of the point m itself, then the above quantities are given functions of the 
variables and their differential coefficients. 

The most interesting case of the latter kind is that in which the motion of the new 
axes is assumed to satisfy the equations 



CO, 



^ = -S, . . (^.) 

X y z 

which express the condition that the instantaneous axis of rotation (of the moving 

axes) always coincides with the radius vector of the moving point m *. 

The radius vector traces, in fixed space, a certain conical surface. It also traces, 
with reference to the moving axes, another conical surface ; and we might always 
assume as one of the conditions defining their motion, that this latter should be any 
proposed surface ; that is, we might assume that the new coordinates x, y, z should 
always satisfy the equation <p{x,y, z) = Q, p representing any given homogeneous func-» 
tion. If to this last assumption we add the two conditions expressed by the formula 
(m.), we further assume that the conical surface traced by the radius vector with refer- 
ence to the moving axes, rolls upon that traced infixed space. 

Suppose, for example, we assume for the equation <p{x,y,x)=0, simply z=0. This, 
with the conditions (&>.), will express that the radius vector is always in the plane of xy, 
and that this plane rolls upon the conical surface traced by the radius vector in fixed 
space. We may then say that the plane of xy is the " plane of the orbit," and that 
the axes of xy, or any lines fixed with reference to them in their plane, are " fixed in 
the plane of the orbit-f-." 

* See Jacobi's first letter to Professor Hansen (Cbellb's Journal, vol. xlii. p. 21). This letter appears to 
refer to some unpublished (?) results of Professor Hansen, which may possibly be similar to those of this 
article. 

t The student of elementary treatises is, I believe, always left to find out for himself what this expression 
means, or ought to mean. 
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Now since we are supposing the equations to belong to a dynamical problem, we 
have u=mx', v=.my', w=mz', and if we substitute on each side of these equations 
the values derived from (85.), art. 72, we find easily 

relations which are true on all suppositions as to the motion of the axes ; but the 
assumptions (».) reduce them to 

u=mx', v=m.y', w=m.z'. 

The further assumption z=0, which involves z'=0, gives w=0, and also (by the 
equations (».)) »2=0. Thus the equations (86.) are reduced to 



x'=^, 
du 


u' 


= . 


_dZ 

dec 


^ dv 


«': 


— . 


dZ 
dy 


0=f, 

dw 


0; 


= . 


dZ, 



dZ dZ 

where in -— and -=-, z and w are to be put =0 after the differentiation. It is to be 

dw ds '^ 

observed also, that the values of u, v, w above given reduce the three first of the 

equations (86.) to the form u=m ■—, v=m ^, w=m ~, from which it is evident that 

du dv dw 

where U is the original force- function, expressed in terms of the new variables. This 
depends only upon the conditions (&>.); but in the case now considered we have also 
w=0. 

Let the origin of coordinates be the Sun ; m a planet disturbed by another planet 
m, whose original coordinates are x^, y^ z„ and new coordinates x^, y^, z, ; also let 
x^+y^+z'=r^ x'^f-^z^=r', xf+yf+zf=rf, a;,^+3/f+2^=rf, and (x,-x)^+(y^-y)^ 
+ (z,— z)^=^^; we have evidently r^=r^, rf=r^, and 

(^-x)^+(y-y)'+(z-z)^=(^-*)^+(2/-2/)^+(:2-^)^ 

and xx,+yy,4-zz,=^x,+3/«/^+s;2,. 

We have then originally 



U= 



mu. . f\ xx^ + yy, + zz,\ 



where («-=m+mass of Sun ; and it is evident that this expression preserves the same 
form when expressed in terms of the ne^v coordinates, and also (which is essential to 

the validity of what follows), that ^, &c. are the same whether the differentiation 

2x2 
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be performed before or after the substitution for x, y^, z,, in terras of a?„_^„ ;s,. If then 
we put z=0 and 

R=«,(l-=±n.), 

we shall have, from the fourth and fifth of the transformed differential equations*, 



^ r^^ dy 



(87.) 



from which it is evident, that, assuming the motion, of m, to be known relatively to the 
new axes, the variations of the four elements of the orbit of m which determine the 
dimensions of the orbit, its position relatively to Imes^xed in its own plane, and the 
time of perihelion passage, will be expressed in terms of the differential coefficients 
of R in the same way as if the plane of the orbit were fixed. But the motion of the 
node of the orbit upon the fixed plane of xy, and its inclination to that plane, must be 
determined by means of the last of the differential equations, as follows : that equa- 
tion gives 

or if we put — O for the term multiplied by ?wm, in the value of U given above, 

and z is to be put =0 after the differentiation, which reduces the above to 

/'dD,\ 

Let i!Wo+»f =«^ so that a is the angular velocity of the plane of the orbit about the 
radius vector; then (observing that M=m«', &c.) we have 

Wq fljj a voOq — MtOj 
X y r m{xy' —s^y)' 



whence va>a—uoDi=-y(xi/' — afy) 



r /d!Li\ 

and therefore m«=: — ™=_(^^j 



r /d£V\ 

^/<j.a{\-(Fy 

which gives a, in terms of the four elements referred to above, and of i. And if we 
put < for the inclination of the orbit to the plane of xy, v for the longitude of the node 



* These equations (87.) have been obtained in a diiferent way by Mr. Bkonwin. Camb. Math. Journ. 
vol. iv. p. 245. 

t See below, art. 80. 
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referred to the axis of x, and (3 for the angle between the axis of x and the node, the 
usual formulse of rotation give 

/=&)„ cos (3— »i sin |3 ■ 

c' sin <=(yo sill iS+cwi cos/3 [• (88.) 

|3'=cw3— i''cos< 

If in these expressions we put 01^=- a,, a>i=^oi,, 0)^=0, and call ^ the angle between 

the radius vector and the node, so that <WoCos(3— »iSin(3=rcos&, a„sin|3+»jC0sf3 
=rsin&, we obtain finally 

;'=acos3^, v'=:oi—. — , 3'= — a sin & cot/, 



sint 



in which the expression above given for a is to be substituted. 

The actual value of (t-J is?wm,«, r-j—ggj, which (since z,, r,, &c. are supposed 

given in terms of f) may be expressed in terms of the four elements first mentioned, 
and t. 

I propose to consider the transformation of the differential equations of the planetary 
theory in a more general manner in the following section. At present I shall add 
some remarks on normal transformations in general. 

74. Theorem. If 

^^=1^1 y^=-r., ........... (89.) 

be a system of 2w simultaneous differential equations, where 

Z=/(.r„ i/„ x^,y^, ...p, q, r, ..., t), 
and ^9,9,^, ... are also explicit functions of«i,&c.,3/i,&c. and t, but are exemptfrom differ- 
eniiaiion in taking the differential coefficients -r-' -j-' ^^^ '^ these equations be trans- 

formed by a normal substitution of new variables li, &c., rji, &c. (art. 62, equation 
(73.)), then the transformed equations are, as in art. 64, 

., d-<V , dT, , d-^ dZ 
drti drii d^i rff,- 

in which Z is expressed in terms of li, &c., >?i, &c., but the differentiations with respect 
to |„ fii are performed before the substitution of these variables in jo, g, r, &c.; in other 
words, p, q,r, ... are still to be exempt from differentiation in forming the differential 
equations. 

This may be proved simply by repeating the reasoning of art. 64. The only dif- 
ference is, that in the term 

^ f c?Z dxj , dZ di/j\ 
'[dsCj drii dyj dr)ij 

dZ dZ 
the differential coefficients ^, ^ are now taken only so far as Z contains Xpi/j inde- 



328 PROFESSOR DONKIN ON THE 

pendently oip, q, r, &c. ; and therefore the term represents the diiFerential coefficient 
J- taken so far as Z contains % independently of p, q, r, &c. The same reasoning- 
applies to the corresponding term in the value of ;j-. The theorem is thus established. 

It is evident that it may be combined with that given in art. Jl, where other func- 
tions analogous top, q, r, ... are introduced by the modulus of transformation P. 

If we call the form of the system of differential equations (89.) canonical when the 
differentiations of Z with respect to Xj, &c., y^, &c. are total, we might call it pseudo- 
canonical when Z contains functions of x^, &c., i/i, which are exempt from differentia- 
tion in forming the differential equations. 

In like manner, if we call a transformation of variables normal, when the differen- 
tiations of the modulus P (equations (73.), art. 62.) with respect to li, &c., y„ &c. are 
total (as in art. 62.), we might call the transformation pseudo-normal when P contains 
functions of the variables which are exempt from differentiation in forming the equa- 
tions of transformation (as in art. 71-)- 

Adopting these designations, we may enunciate the following general theorem of 

transformation : — 

Theorem X. — If a pseudo-canonical system be transformed by a normal or pseudo- 
normal substitution, the transformed equations are also pseudo-canonical, and may be 
formed by the rules applying to normal transformations of canonical systems, provided 
that the functions which are originally exempt from differentiation with respect to 
the variables, be continued exempt to the end of the process ; but if such functions 
occur in the modulus of transformation P, they are subject to total differentiation 

with respect to t \n iovimng the ierm -jr- (See art. 71-) 

[With respect to this theorem there is one important remark to be made. If u, v 
be any two functions of a?,, &c.,i/„ &c. (with or without p, q, r, ... and t), the equation 

/du dv du dv\ /du dv du dv\ 

\dyi dxi dxtdyj~ \dr]i d^i d^i drij ^ ' '^ 

is now only true on condition that the substitution of the actual values of jt?, q,r, ... 
in terms of the variables be not performed till after all the differentiations.] 

75. The theory of the variation of elements affords an interesting example of the 
theorem given in the last article. Consider the following system of differential 

equations, 

,_dZ dD, , dZ dD, 
'^'— %+^' 2/.= -^-^' • • • (90.) 

where in ^> j- the differentiations are total, but H is supposed to contain functions 

of Xj, &c., t/„ &c., which are exempt from differentiation in forming the above equa- 
tions. The system 

,_dZ ,_ dZ 
' dy/ y* dxi 
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is canonical. Let us assume then that a complete set of normal integrals a„ ... a„, 
b^,... b„ of this latter system is known, so that we have 

ai=<Pi{Xi, &c., y„ &c., t), hi—xlx^, &c., y^, &c., t). 

The assumption of these last equations to represent the solution of the complete 
system (90.) is simply a transformation of variables, {a^, &c., b^, &c. being- the new 
variables) ; it is also a normal transformation, since the equations connecting the 
new^ and old variables may be put in the form (see Theorem VII. art. 49, and art. 62.) 

where X (the modulus of transformation) is a function of jPi, ... x^, a^, ... a„, t. The 
function "*" of art. 62 is now obtained by expressing — -^ in terms of «„ ...&„... ; 

but since Z is — -tt expressed in terms of a?„ &c., y^, &c., it follows that when Z is 

expressed in terms of the new variables «„ &c., b^, &c., it becomes identical with W. 
Now if the process of art. 63 be followed, mutatis mutandis, it will be seen that in 
the present case we obtain 

^i— db~ db~^\dyj dbj' 

in which expression the first two terms destroy one another, and the remaining term 
is evidently the differential coefficient of Q with respect to bi, taken so far as Q con- 
tains bi independently of those functions which were exempt from differentiation in 
forming the original differential equations (90.). Similar reasoning applies to the 
expression for &■. 

As this result will be useful, I shall enunciate it separately as 
Theorem XI. — If the original system of differential equations be formed by treating 
certain functions, ;7, q, r, ..., contained in the disturbing function Q., as exempt from 
differentiation with respect to x„ &c., y^, &c., the equations which determine the 
variations of any set of normal elements a^, &c., b^, &c. 

' dbi ' dtti 






on condition that j», q, r, ... be treated, in forming these equations, as exempt from 
differentiation with respect to a^, &c., b^, &c. 

[It is important to recollect, that after these equations are formed, p, q, r, &c. are 
to be expressed in terms of a^, &c., by, &c., and in the integration of the system 
a„ &c., &„ &c. are to be treated indiscriminately as variables, whether they originally 
entered through p, q, r, ... or not]. 

The Theorem XI. may also be immediately obtained from the general equations 
(E.) of art. 52 (in which it is to be remembered that Z includes the disturbing 
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function). The above method of deducing it is given as an additional illustration of 
the general theory of transformation. 

If, instead of the normal elements a^, &c., b^, &c., we employ any other elements, 
Ci, C2, &c., vrhich can be expressed as functions of the former, the formula (69.) of 
art. 54 w\\\ still be applicable, with the condition that p, q, r, &c. are exempt from 
differentiation with respect to c^, c^, &c. 

Section VII. — On the Differential Equations of the Planetary Theory. 

76. The differential equations which determine the motions of a system of mutu- 
ally attracting material points relatively to one of them, do not, as is well known, 
naturally present themselves under the canonical form (I.), art. 49. It is possible 
indeed to reduce them, by different artifices, to that form; but it seems doubtful 
whether any practical advantage is gained by doing so. When the ordinary method 
is followed in the case of a planetary system referred to the sun, there is a distinct 
disturbing function for each planet ; but it is easily seen that the usual expressions 
for the variations of the elements hold good, not merely for each planet on the hypo- 
thesis that the motions of the rest are known, but as a complete and rigorous set of 
simultaneous differential equations involving all the elements of all the orbits, and 
their differential coefficients with respect to t (and containing of course also t ex- 
plicitly). It does not appear that we are practically farther from the attainment of 
the rigorous integration of this system, than we should be if it had the canonical form, 
as it might be made to have if it were derived from an original system of that form ; 
and so far as the development of the disturbing functions is concerned, the most 
troublesome part of them, which is that depending on the mutual distances of the 
planets, is not likely to be got rid of by any conceivable artifice. 

However this may be, all that I propose to do in the present section is to take the 
original differential equations in their ordinary form referred to rectangular axes 
passing through the sun and parallel to fixed directions, and then to exhibit in a 
general manner the effect of a transformation to new rectangular axes, still passing 
through the sun, but changing their directions in space according to any arbitrary 
law. 

77. Let M be the mass of the sun, and m, m^, m^, &c., the masses of the planets ; 
and put M-\-m=(M, M+m,=;M,p M-f m,,=|!A^,, ... M+m^o=iJt>(i). And, referred to the 
original axes, let x, y, z be the coordinates of m, x,, y„ z, of m,, ..., x^^^, y(;j, Z(y of m^^^. 
Also let R, Rp ... R(,) have their usual significations, so that 

^—•^l '^W ^ OT(.-) (xx(i) + yy^i) + zz(i)) 

i((x,)-x)«+(y(O-y)^+(z«-z)0* (^h+n^ + ^fi))^ ' 

the summation extending to all the planets except m. 

Then if we put xf,)-|-yf,)-l-z|,=rfi„ the original differential equations of the second 

order are such as x'li^+ijj^=-^. Let u,.,, Vf^j, W(i) be the variables conjugate to 
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X(j„ y(i„ Z((). Then instead of the original system of differential equations of the second 
order, we have the following system of the first order : 



X(.)_.__-, y(.) 



rfU(„. ''» dw^{ --« ^w« 



in Which Z,,=2-^(uf,+vf,+wf,)-»«c.-)(^+R«)- 

These equations are not of the canonical form, because R(j) is not the same for all the 
planets. But it is easy to put them in the pseudo-canonical form (art. 74.), a process 
which is not necessary, but saves trouble by bringing them under the operation of the 
general rules of transformation established in former articles. 
In fact, if we take 



\ 2ot«) r(i) / 






:)X»)+y(i)y(j) + zwz:o-) , X(i)X(j-)+y(,-)y(^-) + Z(i)Z(y) 



1 



.)■ 



( (X(.-) -xy))^ + (yra - yy))^ + (z« - zy))^)*; ' 

where the summation in the first term extends to all the planets, and in the second 
to all their binary combinations, and a horizontal line placed over any letter indicates 
exemption from differentiation, we shall have 

with similar equations for Y\i>i Z(j), V(,), w'^j,. 

[The terms -^ Sec. are only written for the sake of uniformity, being really =0, 

since Q does not involve U(,), V;;), w^i^.'] 

In these equations Z and Q are the same for the whole system, and the differentia- 
tions of Z are total; but those of Q are restricted to the quantities not marked by the 

horizontal line, so that ^7 is really the same thing as —^(v,-^- 

78. Let us now refer the whole system to new (moving) rectangular axes, whose 
position at any instant with respect to the original axes is defined, as in art. 72, by 
the variable direction-cosines X„, &c. Let x^i^, y^^, z^-, be the new coordinates, and 
«(()> %j ^(i) the new variables conjugate to them. The transformation will be effected, 
as in art. 72, by taking for the modulus 

P== 2((x,X(i) -t- A,2/(i)-H X22(i,)U(j, 4- (po%) H-pi3/(i) +P2«w) V(i) + ( V(i) + ^3/(0 + <'22«) W(i,), 
and the result will be as follows : put 

a=Q-f-«o2(Z(,-,y(j,—3/«,W(.))-|-5,2(%W(j)-S(..)M(i))-|-i82(3^(ott(i,— %,?;(„) 

MDCCCLV. 2 Y 
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(where Q is expressed in terms of the new variables, and a^) '^u '"a (the angular velo- 
cities of the moving system of axes about the three moving axes themselves) are 
marked with the horizontal line to show that these quantities are exempt from diife- 
rentiation in forming the following system of differential equations, though they may 
be functions of the variables) ; then the system (91.) is transformed into 

with similar equations for ?,'•, &c. 

In these equations Z is to be expressed in terms of the new variables ; and it is 
evident from the original form of Z and Q, that when so expressed, these two quan- 
tities are the same functions of the new variables that they were of the old, and in- 
volve (see art. 74.) the quantities exempt from differentiation in the same way*. 
Thus the transformed system (92.) contains no terms explicitly depending upon the 
motion of the axes, except those introduced by the three terms multiplied by u^, a^, % 
in the value of O given above; and the addition of these terms constitutes the only 
difference between the form of the old and of the new system. 

79. We may now apply the method of the variation of elements to the system (92.) 
as follows : — 

The system obtained by omitting the disturbing function Cl, namely, 



, , dZ ^ . AZ ' , dZ ^ 

%)+^=0' ^»+^y«=0' «^«+^,=^ 



(93.) 



is canonical, and consists simply of the aggregate of the equations representing, for 
each planet, undisturbed elliptic motion about the sun -f~ (relatively to the new axes 
of coordinates). 

The integrals of these equations may therefore be expressed in any of the usual 
forms. We will suppose that the elements chosen are 

a,e, ■m, (s), i,v, 

with significations corresponding to those given to the same symbols in art. 55. 
These letters unaccented will apply to the planet m, and a„ e^, &c., a^^, e,,, &c., «(j„ e^iJ, 
&c., to the planets m,, to,,, m^^^, &c. 

The definitions of the elements a, e, ts, &c. are their expressions in terras of the six 

* Since the direction cosines X^, &c. are exempt from diiferentiation in forming the equations connecting the 
old and new variables from the modulus P, they continue exempt throughout. (Theorem X. art. 74.) Hence 

and similarly for the rest. 

t Z=:Sl H — j, the summation extending to all the planets. 
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variables x, y, z, u, v, w*, and t, and the same expressions continue to be their defi- 
nitions when they become the variable elements of the disturbed motion. 

Now the general formulse for the variations of the set of elements here chosen 
have already been given in art. 55 ; for it is evident that the process in the present 
case would merely be a repetition, for each planet, of the process there employed. 
Here however we are to use in every case the disturbing function fi given in art. 78 ; 
but if we observe the effect of the marks of exemption, it will be evident that, for the 
planet m, the only effective part of O is 

a^{zv~yw)-{-aii{xw—%u)-\-&)2{yu—xv)—mR; 

and similarly the effective part of Cl for any other planet will be given by suffixing 
the corresponding number of accents to x, y, z, u, v, w, m, R. 

Now if we put A„, A„ A^ for the terms multiplied respectively by ao> '^n % ^^ the 
above expression, we have, by the definitions of the elements, 

{Al+A', + Alf=mViMa{l-e') 



A2=— ?n\/j«-a(I — e^). cos i, Ax=:m\/ ^a{\ — e^) . cos v sin / 
A„= —m\/^a{ 1 —^) . sin v sin <, 
so that the disturbing function, so far as m is concerned, becomes 

— m{R-l-\/ja/«(l — e^).(ft'oSini'sin<— •ftiiCosi'sin<+*'2COS<)} (94.) 

Consequently, since the expressions in art. 55 were obtained by taking — wR for the 
disturbing function, we have merely to add to them the additional terms derived 
from the part added to R in (94.). Performing the differentiations, and omitting 
afterwards the symbols of exemption over'<yo5 ^\, <>>% which cease to be of any use, we 

obtain, after obvious reductions, the following simple results : if -^, -^i &c. represent 

those parts of the differential coefficients of a, e, &c., with respect to t, which depend 
upon the motion of the axes of coordinates, then 

dt ^' dt ^ 

"hU) Bar . </ . x 

-^ = — -- = tan 2(<*'i cos v—co^ sin v) —u^ 

-Lz= — {u^ COS v-\-coi sin V) 

-l=:—co.2-\- cot/ (»(,sin)'~», cosv) 
where it is evident that we may write s instead of (s) (see art. 55.). 

* Not their expressions in terms of ce, y, z, x' , y', z' ; for though these are equivalent in the undisturbed equa- 
tions, they are not here equivalent in the disturbed equations, and therefore the general theory, which assumes 
the former mode of expression, is not here applicable to the latter. 

2 Y 2 



(95.) 
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The complete variations of the elements are then found by adding to the terms 
just written the expressions given in art. 25. 

It is easily seen that the expressions (95.) might have been deduced from geome- 
trical considerations alone, if we had been at liberty to assume beforehand that the 
mechanical and geometrical parts of the variations might be calculated separately; 
the former as if the axes were at rest, and the latter as if there were no disturbing 
forces. It would not, I believe, be difficult to establish by ct priori and simple 
reasoning the validity of such an assumption, and then the above results would only 
serve as a verification of the method which has been employed to obtain them. 

80. In order however that no obscurity may rest upon the interpretation of the 
formulae obtained in the last articles, it is necessary to consider the physical (or rather 
geometrical) meaning of the elements a, e, &c., which we have so far only defined by 
means of their expressions in terms of the variables x, y, z, u, &c., and to ascertain 
what relation they bear to the elements similarly defined by means of the original 
variables x, y, &c., which refer to axes whose directions are invariable. 

The relations between the variables ((85.), art. 72.) give immediately 

*'+/+ 2;'=x^+y'+ z' 

M^+W^+M)^=U'+V'+W' 

ux-\-vt/-^zw:=ux-\-vy-{-zw; 
and if we put yw—zv=^A, zu—xw=B, xv—yuz=C 

yw— zv=A, zu— xw— B, xv— yu=C, 
we find, by virtue of the relations piVj— yia.2=Xo, &c., the following equations : 

^zriX^A + ^oB + VoC 

5=X,A+|«-,B+j',C 
C=X,A+/«,,B+v,C 
and ^^+5HC^=A*+B^+e. 

Now A(=yw— zv=m(yz'— zy')) is the projection on the plane of yz of the areal 
velocity of m (relative to fixed space) multiplied by the mass, and B, C have analogous 
meanings ; hence it is evident from the above equations that A, B, C are the pro- 
jections on the three moving coordinate planes of yz, zx, xy of the absolute areal 
velocity of m relative to fixed space, multiplied by the mass. (The projections of the 
areal velocity relative to the moving ao^e* would heyz'—zy\ &c., which are not pro- 
portional to yw—zv, &c., since u, v, w are not the same as mx', my', mz', except on a 
particular hypothesis as to the motion of the axes. See art. 73.). 

Inasmuch as the definitions of the elements a, e, i involve the variables only in the 
forms x^+y^-^-z"^, M^-fw^+w^ A, B, C, it follows that these three elements are 
respectively the semiaxes, excentricity, and inclination to the plane of xy, of the absolute 
osculating ellipse of the orbit infixed space. Thus the instantaneous ellipse, relatively 
to the moving axes, is of the same dimensions and in the same plane as the true 
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osculating ellipse ; and it only remains to show that it coincides with the latter in 
position, for which purpose we must prove that it touches the true orbit. (It does 
not in general touch the relative orbit, because mx', wy\ mz' are not in general the 
same functions of the elements and t that u, v, w are.) 

81. If we suppose the coordinates x,y, 2 of m expressed in terms of the elements 

a, e, &c. and t, and denote by -^) &c. their differential coefficients taken so far as t 

Git 

appears explicitly in these expressions, then ^> ^3 ^ are proportional to the direc- 
tion cosines, relatively to the axes of x, y, %, of the tangent to the (relative) instan- 
taneous ellipse. And therefore the direction cosines of this same tangent referred to 
the fixed axes of x, y, z, are proportional respectively to 

. dx , ., dii , . dz dw , dy , dz dx , dy , dz 

'-H-^Ht-^'^Ji ^^Tt^^4+^^Jt' '"Tt-^'^t'^'^-di 
On the other hand, the direction cosines of the tangent to the absolute instantaneous 

ellipse, referred to the fixed axes, are proportional to -^3 -J-, -| (the differential 

0/t wt at 

coefficients of x, y, z, taken so far as t appears explicitly in the expressions for those 

variables in terms of the elements of the absolute ellipse). The identity of the two 

tangents will therefore be established, if we can show that 

dx - efe , , dy , . dz o 

Now -^ is that part of x' which does not depend upon the disturbing function ; i. e. 

(equations (92.), art. 78.) we have identically 

dx dZ dy dZ dz dZ . 

dt du dt dv dt dw 

J • i-i dx dZ dy d!L dz d!L 

and, m like manner, t-=t~' -t-==-t-' -r=-T-' 

dt dn dt dv dt dw 

where Z is the same as before, but expressed in terms of x, y, z, u, v, w, instead of 
x,y,z,u,v,w. Now let the latter set of variables be expressed in terms of the former 
by the formulae of art. 72, and we have 

rfZ dZ du dZ dv dZ dw 

du du du' dv du''dw du' '' 

, du dv dw , . s „ 

but 5ii=^»' ^11=^" rf[r=^= ((^^•)' ^^'^' 72-)3 &c., 

I I r ^x dx , ^ dy , dz a 

and therefore lt=^odt+^^li+^Hi' ^^•' 

as was to be proved. 

82. It follows, then, that the mode of treating the problem adopted in the pre- 
ceding articles is equivalent to representing the motion of each planet by means of 
the true osculating ellipse of its actual orbit (relatively to the sun) in fixed space. 
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The definitions of all the elements (relative to the moving axes) in terms of the six 
new variables x, y, z, u, v, w, have the same form as those of the corresponding 
elements (relative to the fixed axes) in terms of x, y, z, h, v, w. The two relative 
elements a, e are the same as the corresponding absolute elements ; / is the inclina- 
tion of the plane of the ellipse to the moving plane of xy, and i> the longitude of the 
node reckoned from the axis of a;; and since the place of the body in the ellipse is 
evidently the same, the relations between the remaining elements srand (g) (or s) and 
the corresponding absolute elements are purely geometrical. 

Comparing these results with those of art. 79, we see that the independence of the 
formulae for the mechanical and geometrical variations of the elements of the true 
osculating ellipse is completely established. 

83. In all that precedes, the three variables &)a, a)i, tu^ (the angular velocities of the 
system of moving axes about the axes themselves) are entirely arbitrary ; they may 
be either explicit functions of #, involving only determinate constants, or they may 
depend in any way upon the relative or absolute elements of the orbits of any or all 
the planets, and their differential coefficients with respect to t. In the case in which 
the expressions for »„, »„6;2 involve only the relative elements, when these expressions 
are introduced in the formulae (95.), art. 79, and these formulae completed by the 
addition of the terms in art. 55, and when the corresponding sets of equations are 
formed for each planet, we obtain a set of simultaneous differential equations in- 
volving all the elements of all the orbits, and their differential coefficients with 
respect to t. The integration of these equations would determine all the elements 
as functions of #, and thus the motion of all the planets, relatively to the axes of 
coordinates, would be known. Lastly, the motion of the whole system, relatively to 
fixed space, would be found by integrating the system of equations 

L'=&'oCosX— WiSinX 1 

N'sinL=^oSinX+^jCosX ^ • • • (^^•) 

X'=*>,-N'cosL J 

where u„, Oi, w^ are now given functions of t, and L is the inclination of the plane of 
xy to that of xy, N is the longitude of the ascending * node of the plane of xy, 
reckoned from the axis of x, and X is the longitude of the axis of x, reckoned upon 
the plane o{ ry,from, the node, in the direction of positive rotation. 

In the case in which a^, u^, u^ cannot be expressed in terms of the relative elements, 
the integrations which determine the relative motion of the system cannot be sepa- 
rated from those which determine the position of the axes in fixed space ; but the 
equations (96.) must be considered simultaneously with the other differential equa- 
tions of the problem. 

* Ascending relatively to a positive rotation, i. e. from x to y. 
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Application to the Problem of Three Bodies. 

84. I propose to exemplify the preceding principles by applying them to the trans- 
formation of the equations which determine the motion of three mutually attracting 
bodies, considered as material points. Let them be called the sun and two planets, 
and let the origin of coordinates be placed at the sun, and the notation be the same 
as in arts. 77-79, so that M, m, m^, are the three masses, M.-\-m=^, M4-?«j=i"'„ &c. 
Let the elements be chosen as in art. 55, and longitudes be measured, as before, 
along the plane oixy (from the axis of*) as far as the node, and then along the plane 
of the orbit. 

Putting for convenience mR=n, w,R,=np we have 

VL-=.mm,{{r^-\-rf—2rriCO%y)~^^—rr~^co^X} 
ll,=WOT^{(r^+''/— 2rr,cos y)--^—r^r~^ co^ %}, 

where % is the angle between the two radii vectores. Let 1 be the mutual inclina- 
tion of the planes of the two orbits, and let the angular distances, on the planes of 
the orbits, between their ascending nodes on the plane of xy and their line of inter- 
section, be respectively v, v, ; so that 6—v—i>, d^—v—v^ are the angular distances of 
the two radii vectores from the line of intersection (which we may call, simply, the 
line of nodes). We shall then have 

cos%=cos(^— u— p)cos(^^— y^— !',)-l-sin(^— y— i')sin(^j— y,— i';)cosI| /„^ ^ 

cosI=cos/cos/,+ sin<sin;;Cos(v,— v)* i 

and u, tf, are functions of /, ^^, v^—v, determined by the equations 

cotysin(V(— v) = — cot;jSin/-|-cos(v,— i')cos<"l .j,„ ^ 

cotyjSin(i',— i') = cot/sin/,+ cos(v^— >')cos<^ J 

Now considering Cl as expressed, on the one hand, in terms of r, r,, 6, 6^, v, v^, i, i^, ami 
on the other in terms of all the elements and t, we have rigorously, as may easily be 
proved in the usual way-j~, 

di ds ~^dm' di) (/=,' dia^ 

* The arrangement referred to here and in the folio wing 
articles will be made clear by the accompanying diagram. 

f Since the values of r and fl in terms of the elements and 
i are of the forms 

r—f{fnit-\-i—'^'), i=fttdt+»+<p(/ndt+t—'m), 

, dr dr ^ d$ , di ^ 

we have + =0, —+—=1; 

as am as aw 

from which the equations in the text follow immediately. It 
may be as well to remind the reader who may happen to 
recollect the note to the Astronomer Royal's tract on the 
Planetary Theory (p. 91, ed. 1831), that that note refers to 
a different way of measuring longitudes. 
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with similar equations for H^. In what follows, I shall, as an abridgment, employ 

de diss 



3 J 

the symbol E to denote the operation j+;7-> and in like manner I shall put E, for 



A+A, so that 



Ea=^+^,&c (99.) 

ds dvT ^ ■' 



Since r and 6, when expressed in terms of the elements and t, do not contain /, /„ v, »,, 
we have 

dH dCl dcosx dD, dD, dcos^ 

dv dcosx dv di dcosx <^' 

. ^cosx_ { ^^du\ d cosx du, dcosx 

+ sin / sin i, sin (6 — v—v) sin (^,— y,— v,) sin (v—v), 

in which expression the values of j) -^ are to be obtained by differentiating the equa- 
tions (98.). In like manner we find 

dcosx '^" dcosx__dtJi dcosx 

di di di di fdi/ 

+ sin {d—v—i>) sin (^/— «;— 0(— sin/ cos /,+ cos i sin i, cos (>',—»')). 

Analogous expressions may be found for -j-', -—'• 

85. Hitherto we have assumed nothing concerning the motion of the axes of 
coordinates. Let us now however take as a first assumption that the plane of xy shall 
always pass through the line of nodes. This implies the conditions* 

"/="> w=0, W/=0, 
and consequently I=<^— /. 

Introducing these conditions in the expressions at the end of the last article, and 

* The legitimacy of the following processes will be apparent to the reader who shall have followed the 
general reasoning of former articles, though probably not to others. In either case it may be useful here 
briefly to recapitulate the principles now to be applied. The results of art. 79, and in particular the expressions 
(94.), (95.), were established independently of any assumption as to the values of ctf^, Wj, ui^, which are perfectly 
arbitrary. We are therefore at liberty to assume that W(,, u^, w^ are such functions of t that any three func- 
tions of the variables shall constantly =0, Thus the first assumption made in the text is that y^— y shall con- 
stantly =0. But since such assumptions may cause (as in this case) certain elements to disappear from the 
expression of CI, it is necessary to perform the differentiations of ii with respect to such elements /rs# ; the dif- 
ferential coefficients may or may not vanish, on afterwards introducing the assumptions ; if not, we find ex- 
pressions for them in terms of differential coefficients with respect to other elements which do not disappear ; so 
that instead of differentiations which ought to be performed before the assumed conditions are introduced, we 
have finally only such as may be performed afterwards. The expression " disappear" does not, it must be 
observed, necessarily mean the same as "vanish." Thus the condition y^— »'=0, causes v and V; both to dis- 
appear from O, but does not of course imply that they both vanish. 
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dQ, 



observing that -^=EQ, &c. (see equation (99.)), we have evidently 



dD, 



di' 



U-(,+(|))E.-(|)E,a 
.(S)E.-(f)Ep-f; 

and it only remains to find the values of \j^), &c., namely, the values of ^j &c., 

which correspond to the assumed conditions v^—v=Q, u=u^=0. 

Now in any spherical triangle of which a, b, c are the sides, and A, B, C the oppo- 
site angles, if a be considered as a function of c. A, B, by virtue of the equation 

cot a sin c= cot A sin B + cos c cos B, 
we have by differentiation 

■T-= (sin a)^(cot a cot c+ cos B) 

= cos a sin « cot ^+(1 — (cos a)^) cos B= cos B-|- cos a cot C sin B 

da /sinaX^sinB sinasinB 



da /sma\^ 

SS~"\sinAy 



sine 



sinC 



^ = _ (!i!Lf)! cot A cos B + cot c(sin af sin B = sin a cot C ; 
dB sine 

and if the sides of the triangle be indefinitely diminished, these expressions become, 

in the limit, 

da sin(B + C) da ^ da^ ^ 

Tc sinC ' dK ' dB~ ' 

provided the angle C do not vanish. 

If these results be applied to the triangle of which the sides are v, v,, v^ — v, and the 

opposite angles 5r— /,, <, I, it is easily seen that the values of [lu, &c. are as follows : 

(du\ sin«( (^^\ ^^°'/ (^^i\ ®'^'^' i^^\ ^''^' 

\dv) sinl' \dvj sini' \dv / sini' \dvj sin I 

provided I be not =0. We have therefore 



dv 



■^l^(sin/^.EQ+sin/.Ep) 



m 
dr 



da 

-di' 



(100.) 



exactly in the same way we find 



dVf I ' sm 



j(sin/.E,Q^+sin<,.EQ^) 



dUf dD,i 
dt, dl ' 



(101.) 



MDCCCLV. 



2 Z 



(102.) 
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The variations of the elements will now be found by introducing the above values of 

da da, dP, dP,, 
dv dvi dt dii 

in the expressions given in art. 55, and completing these expressions by the addition 
of the terras (95.), art. 79. But the angular velocities aio, a^, u^ are no longer wholly 
arbitrary, since we have made one assumption concerning the motion of the axes, 
which implies one relation between these quantities and the elements. In order to 
determine u^, m^, u^ completely, it will be necessary to make two more assumptions ; 
but first we will investigate the relation already implied. 

86. The complete expressions for v, v^, obtained in the way mentioned in the last 
article, from arts. 55 and 79, may be written in the following form : put for brevity 

and put |tA*a-i for n, and i/^^afi for n^ ; then 

<=^, ^'+cotv(«oSin.,-«,cos.,)-<., 

and if the latter equation be subtracted from the former, and the conditions 

_ m__dD, dD,i_dQ,i _T 

be introduced, the result is easily found to be 

, . , . _ sini dO,, , sin I, dD> 

{a,,smv-a>,eosv)smI=-j^ -^+~j^^. ..... (103.) 

.This is the relation between a/o, »i and the elements and t, implied by the one 
assumed condition that the plane of xt/ passes through the line of nodes. The angular 
velocity ai^ of the system of axes about the axis of z, is so far left, as it evidently 
ought to be, perfectly arbitrary. 

Q and Q, are now functions of the following elements* only: — 

a, e, (s), m, a,, e,, (g^), w„ I, c. 
And we now have 

cos%= COS (^— J*) cos {&i—p)+8m (d—v) sin (^,— i*) cos I. 

87. The complete ex()ression for /', derived from arts. 55 and 79, is easily put in 
the form 

and, on introducing the assumptions of the preceding articles, this will be found to 
become, after simple reductions, 

p sin 1.;'= — (cos I.EQ+E^Q) — (aocos v+aiisini') jB sin I ; 

* On the difference between (s) and e see above, art, 55. "We cannot strictly call ii a function of « and s. 
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similarly, we find 

jo,sin I./=cos I.EPj+EQ;— (iWoCos c+^i sin f) />; sin I ; 
and since 1'=/—/, we obtain from these 

j9p,sinI.r=_p(cosI.Ep,+EQ^)+;j,(cosI.EQ + Ep), . , . (104.) 
an equation which may be transformed as follows : — since 

W=0> |=0(a«.79), 

we find from art. 55, P'~'d — l"J~ — ^^' 

and similarly jB|=E^Qp whence it is easily seen that the above equation may be written 

in the form 

{pp,cosl)'+pEQ,+p,E,Q=0. ....... (105.) 

If we investigate, from the above expressions for / and /, the values of (p sin /)' and 
(/>;Sin/j)', we find 

sinl.(jjsin;)' = — (cos;,. E+cos/.E,)Q—^ cos; sin !(»„ cos c+sUi sine), 
■ sinl.(/),sin;,)'=(cos;.E,+cos;y.E)Q,— jp,cos;,sin I((y„cosj'+iyiSini'), 

and, adding these equations, 

sin I.(p sin i-{-p,sin ;,)'=(cos ;,.E4-cos ;.E,)(Q,— Q) 

— (pcos;+j9;Cos;,)sinI.(«UoCOSi'+®iSini'). . . . (106.) 

With respect to this equation and (103.), it may be observed that ^oCosv+iWiSin v 
is evidently the angular velocity of the plane of ocy about the line of nodes, and 
iiijCosy— ft'oSin!' is the angular velocity of the same plane about a line in itself per- 
pendicular to the line of nodes ; ori which comes to the same thing, the angular 
velocity of the line of nodes itself in fixed space, estimated perpendicularly to the plane 
ofxy. 

88. The position of the plane o^ xy at any instant has been so far left arbitrary, ex- 
cept that it has been subjected to the condition of passing through the line of nodes. 

As a. further assumption, that which most naturally presents itself is, that the plane 
ofxy should always coincide with the principal plane. By the principal plane I mean, 
of course, that on which at any instant the sura of the projections of the areal velo- 
cities (multiplied by the masses) of the two planets about the sun, is a maximum ; it 
evidently always passes through the line of nodes, and would be the invariable plane 
if the disturbing functions vanished. 

To determine the position of this plane we have (see art. 80.) to express the con- 
dition that ; and ;, are always so taken, subject to the equation ;,— ;=I, that the ex- 
pression mv^fia(F^^^^. cos; +?ra,\/|M<^a,(l—e,^). cos;, shall be a maximum. We will 
put ff for the value of this expression, so that using p and p, with the same nieaning 
as before (see art. 86.), we have 

ff=j9 cos ;-fp,cos ;p 
2 z 2 
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and the required condition of a maximum will evidently be 

j9sin;+/';S"i'/=0*; (107-) 

adding the squares of these expressions, we obtain 

a^=p^+p^+2pp,cosJ, (JOS.) 

which determines the actual value of (t; moreover we have 



sin* siniy sini _ Ciaq \ 

p, p 0- ' ■ • ■ ^ ■'' 

and it is easy to find ffcosi=p-^PiCosI 1 (110.) 

(T cos ti=pi-\-p cos I J 

so that a, sin I, sin/,, cos<, cos;, are all simply expressible in terms of j9, p, and I. 
The variation of c is easily found by means of the equation (105.), which gives 

<r<^=(p,-E,-pE){Q,-Q) (11 J.) 

The equation (103.) now gives (see (109.)) 

o'(<w,cos)'— <yoSin)')=^(Q;— Q) ; (112.) 

and from (106.) we obtain 

<r^sinI(<yoCosc+aiSini') = {(j»+p,cosI)E^+(jo^+j9COsI)E}(0;— Q). . (113.) 

The two last equations determine the motion of the principal plane in space, irre- 
spectively of any arbitrary sliding motion which we may attribute to it in its own 
plane. For they give the angular velocities with which it is at any instant moving 
about two lines at right angles to one another in its own plane (see the end of art. 87.). 
They may be put in another form as follows : — the actual value of H, — O is 



mml~—^jcosx, 



where cosp(^ has the value given above (end of art. 86.); and if the operations indi- 
cated be actually performed, observing that Er=0, E,r=0, &c., and that 

Ecos;^=— ^, &c., 
the results will be found to be 

(7* T \ 

^- — jjsin(^— j')sin(0, — v) 

<r'(<wocos)'+<yisin)')=mm,sinl.(p— ~| x (j»cos(^— j')sin(^^— )')+j9,sin(^— v)cos(^,— )/)). 

Here ^— f, (^,— v represent, it will be remembered, the angular distances of the planets 
from the line of nodes. 

We will assume for the present the condition a>2=^, so that the plane of a?j/ may 
have no sliding motion, but roll upon the conical surface to which it is always a tan- 

* Referring to the arrangement supposed in the diagram, it will be seen that < becomes negative in the case 
now considered. 
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gent plane. If, then, we call ca the actual angular velocity of the principal plane 
about its instantaneous axis of rotation, so that w=\/u\-\-m\, and put j for the angle 
between the line of nodes and the instantaneous axis (which is the line in which the 
principal plane is intersected by its consecutive), we shall have 

«yoCosj'+a,sinj'=iiycosj and ca^Q,o^v — Uf^^\Viv=.ca%\x\j; 
and the above equations give 

(T cot J =J0 cot (^ — I') +jo^cot (^, — v) , 
a result which may also be put in the form (see (109.)) 

sin I cotj=sin /^cot (^— y) — sin < cot (^^— k). 

It is very easy to show by spherical trigonometry that this equation signifies that 
the instantaneous axis coincides with the line in which the principal plane is inter- 
sected by the plane of the two radii vectores*. In other words, we have this theorem : 

The principal plane always turns about the line in which it is intersected by the plane 
of the two radii Vectores. 

It follows of course that the principal plane always touches the conical surface 
described in fixed space (relatively to the sun) by the said line. I think it probable 
that most persons would expect, at first sight, that the principal plane would always 
touch the conical surface described by the line of nodes, which, as has been just 
shown, is not the case. It is perhaps worth while to verify this result by inde- 
pendent reasoning. 

89. Let Roman letters refer, as in former articles, to axes of coordinates having any 
arbitrary;^xe6? directions. Then, putting A=m(yz'— zy'), &c. (as in art. 80.), and 
using §, }}, ^ for current coordinates, the equation to the principal plane is 

(A+A,)|+(B+B>+(C+CX=0; (114.) 

and the line in which this plane is cut by its consecutive is determined by combining 
the equation (114.) with that obtained from it by differentiating the coefficients of 
I, jj, ^ with, respect to t; namely, 

(A'+A;)|-f(B'+B'>+(C'+CX=0 (115.) 

Now from the fundamental equations 

we obtain A'=m(yz"-zy")=y-^-z-^=mmXyz,-zy,)(S-'-r-') 

* Let ^p be the angle between the principal plane and the plane of the radii veetores ; the former plane 
divides the angle I into two parts, of which one is ;^ and the other is — i; and we get two spherical triangles 
which have a common side/, "with adjacent angles *j and vf/ in one, — f and v—^ in the other; hence 

cot(9— >')sin_;=cot»J'sini + cosycos( 
cot(S,— y)sin^'=cot rj/ sin (, + cos^' cos I,; 

and if cot ^ be eliminated between these, the result is the equation in the text. ■ 
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(where 5 is the distance between m and m) ; and in like manner 

A'=— wm/yz^— zy,)(S-'— r-«) ; 

we have therefore A'+Aj— mm/yz^— zy,)(r^~'— r~^), 

with similar expressions for B'+B[, C'+C^; so that, when the common factor 
m,m,{v~^—v~^) is omitted, the equation (115.) becomes 

(yZ/-zy/)^+(zx,-xz>+(xy,-yxX=0, 

which is evidently the equation to the plane containing the two radii vectores. Thus 
the theorem in question is verified. 

90. To return from this. digression : the motion of the line of nodes in the principal 
plane will be given by putting »2=0 in either of the equations (102.), and introducing 
the value of a^sinf— li'iCosv from (103.). In this way we find, after slight reductions, 

pp,sinIJ=:p^cosi,^+pcosi-^, 

in which we may substitute for cost, cos<;, the values given by equations (110,); this 
gives 

ffjop; sin I . v' = {p^+pp, cos I)^ + ip''+PP, cos I)^' ; 

or, if we introduce the actual values of -jj-' "jt^' we find 

pp/= — mm^rr iSin (6—») sin (^(— f) X {p cos / . (§"' - r~^) +J9;C0S /,(^"' — r~^) } . 

It is not my purpose however to enter further into details; and I shall conclude 
this subject by briefly examining the consequences of a slightly different assumption 
as to the motion of the axes of coordinates. I shall suppose, namely, that the plane 
of «y still always coincides with the principal plane, but has a sliding motion such 
that the axis of « alwafys coincides with the line of nodes. 

91. The assumption made at the end of the last article implies the condition v=[) ; 
and »2 will no longer be 0, but must be determined by equations (102.) ; either of 
these gives (putting /=/=0, and reducing by means of (103.), (109.), &c.) 

pp,sml.a)^=p^cosii-jj-\-pcos,i-^, (lib.) 

which coincides, as of course it ought, with the expression given for v' on the former 
hypothesis (art. 90.). The difference is that Q, Q, now no longer contain v. 

The values of ^oj ^i are obtained at once by putting ^=0 in the equations (112.), 
(113.) ; and all the conclusions which were derived independently of any supposition 
as to the value of ai^, subsist as before, when modified by putting i'=0. 

We may add one more equation, which is required in forming some of the expres- 
sions for the variations of the elements ; namely, 

tani^-^i^'^^^-v. (117.) 

2 (T+p+piCOsl 
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This is easily obtained from (109.) and (110.) ; and in like manner 

tan^=-^iHLl . 

2 0" +Pi +p cos 1 

I shall now recapitulate the results of the last supposition, so as to exhibit in one 
view the transformed differential equations of the problem of three bodies. It will 
be as well to repeat also the explanation of the symbols. 

92. Signification of the Symbols. 
M, m, m, are the masses of the sun, and of the two planets, and M+m=jM,, 

a and e are the semiaxis and excentricity of the instantaneous ellipse described by 
m about the sun. 

a, and e^ have the same meanings with reference to m,. 

I is the inclination of the plane of the former ellipse to that of the latter. 

0, 6^ are the longitudes of the two planets, measured in the planes of their orbits 
from the common line of nodes. 

r, r, their radii vectores. 

fs, ns^ the longitudes of the perihelia, measured likewise in the planes of the orbits 
from the line of nodes. 

i, g, two elements such that \ ndt-\-z, \ n^dt-\-%^ are the mean longitudes, where n, n^ 

are defined as usual by the equations n^=i/ja~^, n'^i/jpr^. 

Then I ndt-\-i—Gs is the mean anomaly of m, and r, 6 are functions of the mean 

anomaly and mean longitude given by the laws of elliptic motion. The same is true 
for m„ mutatis mutandis. 
X is the angle between the radii vectores, so that 

cos x= cos ^ cos ^;+ sin ^ sin ^^ cos I. 

Let t be the distance between the planets, so that 

l^=.r^-\-r^^— 2rri cos %. 

Q, Q; are the disturbing functions, defined by the equations 

„ /I rcosy\ ^ /I r, cosvN 

Q=mm.{j ^j, Q-mm,{j-^-^y 

and, when expressed in terms of the elements and t, are functions of 

a, «p e, e^, I, I ndt-\-i, \ ndt-\-i—m, j w^c?^+e^, j w^rf#-l-s, — ra",. 

Jo Jo Jo Jo 

When Q, O^ are considered on the one hand as expressed in this way, and on the 
other, in their original form as functions of r, r„ S, 6,, and I, we have, as applied to 

either of them, 

d d d d d d 

'M~Je~^dw' dS'i~d^i~^d^i 
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p, Pi, a are defined by the equations 

/, /; are the angles between the principal plane and the planes of the two orbits, and 
are given functions of the elements a, «„ e, e^, I, by virtue of the equations 

— jso- sin ;=^,ff sin /,=J9|>, sin I, 
from which we have also 

ffcos/=^+p^cosI, o-cos<,=/?,+pcosI 
Tfor the values of tan --, tan^j see art. 91. j. 

(>)„ is the angular velocity of the principal plane about the line of nodes ; 

jy, the angular velocity of the principal plane about a line in itself perpendicular to 
the line of nodes ; 

% the angular velocity of the line of nodes estimated in the direction of the prin- 
cipal plane. 

Differential Equations of the Problem. 

93. The nine intrinsic elements, as we may perhaps appropriately call them, namely, 

a, a^, e, e,, z, z„ w, w,, I, 

are determined as functions of t by the following system of nine simultaneous differ- 
ential equations of the first order : 






OTjM. Jcos \ dD, \ (Kl^ 
sin I [ p dl 'pi dl J 



'^" ' -da,^ e, ^J V i ej^^ 

»?,/*, Jcos I «?X1, 1 JO 
sin ly Pi dl ~^p dl 

^„^^/-. w«-/l— e^ dD,_ mf/, fcosl ^ i 1^ dD^ 
^ e de sin I [ p dl '^p, dl 

I nfli Vl—e^ dD,i m;/<.,/cosI <?0, , 1 dD,']^ 

ii^inx- — — , _ innl"^ Ti '^p dij 
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PP' «^" I • I-i^.fOS l(js+d^) +-^+^] 

+,{eosl(§+§;)+f+S}- 

94. The only parts of the preceding expressions of which the deduction is not per- 
fectly obvious, are the terms involving I in the values of s', s', w', ot'. They are 
obtained, as has been sufficiently explained, from the expressions in art. 55, to which 

Be 

are to be added the values of ^' &c. ((95.), art. 79.) ; on putting y=0, a;„ disappears 
from the latter; and the values of »„ co^, cos/, cos/^ tan^? tan^ are to be introduced 

(equations (110.), (112.), (116.), (ri7-))- After some rather troublesome reductions, 
the expressions above given will be found. 

In these equations it will be recollected that the mean longitude of m is represented 

by \ ndt+i, and the differentiation with respect to a in ^ is only performed so far 

as a appears explicitly. If we wished that the mean longitude should be expressed 
by nt+s, the only change in the equations would be that the differentiation with 
respect to a must be total ; i. e. must extend to a as contained in n. A similar 
remark applies of course to e,. 

In actual use it would be more convenient to introduce R, R, instead of Q, O, ; the 
latter functions give a rather more symmetrical form to the equations, and are more 
convenient in general investigations. (The relation between them here is merely 
Q=:mR, Q^=OT^R, ; in another part of the paper the symbol Q was used for — wzR 
(art. 55.)). 

95. If the equations of art. 93. were completely integrated, the intrinsic motion 
of the system would be completely determined; that is, we should know at any 
instant the dimensions of the two orbits, the mutual inclination of their planes, the 
position of their axes with respect to the line of nodes, the place of each planet in its 
orbit, and (by (110.)) the inclination of each orbit to the principal plane. 

The position of the system relatively to fixed space would then have to be sepa- 
rately determined as follows : — 

The three quantities »„, &>„ a., (see end of art. 92.), of which the values are ((112.), 
(113.), (116.)) given by the equations 

<r^sinI..„=|(p-hp,cosI)(|^-f^) + (;>,+pcosI)(|+^)}(^-0) 

d 



ff<yi=^(^,— ^) 



.7 sin I .<.,= (^'+ cos l)^-f {^-^ cos l)^, 

MDCCCLV. 3 A 
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would be given functions of t. Then if we call 

J the inclination of the principal plane to an arbitrary fixed plane ; 

Q the longitude of the line of intersection of these two planes, reckoned in the 
fixed plane from a fixed line ; 

N the angle between this line of intersection and the line of nodes; 
we should have (as in art. 83. with a different notation) 

J'rzft/gCosN— »jSin N t 

Q'sin J=:<i»„sinN+»iCosN i (H8.) 

N'=a)2— cot J(&'o sin N+a), cos N)J 
and the integration of this system would give J, Q, N as functions of #, and so deter- 
mine at any instant the position of the principal plane and of the line of nodes, rela- 
tively to fixed space. 

With respect to the motion of the principal plane, the following may be added. 
It has already been shown (art. 88, 89.) that the line about which it is at any instant 
turning, coincides with that in which it is intersected by the plane of the radii vectores ; 
and the values of a>a, wj (see art. 88, putting ^=0 in the expressions there given) may 
be put in the form 

o-X=^w*/ sin 1. 1 4 — 5 1 {p cos 6 sin 6,-\-p, sin 6 cos 6) 
~i — 2) sin ^sin^,. 



ff&i 



If the latter of these be multiplied by cr, and then both sides of each squared, and the 
results added (after putting for c^ on the right its value p'^-{-p^+2pp^cosI), we find, 
observing that cos 6 cos ^^-j- sin 6 sin ^^ cos 1= cos ;;^:, 

(r^\/s»l-\-a)l=mm^ sin I ( ^ — ^j (p^ sin^ S^+p^sin^ ^-\-^PPi sin 6 sin 6, cos %)*, 

an expression which may be further transformed as follows. Let X, X, be the latitudes 
ofrrt,w^(with reference to the principal plane) ; then sin X= sin ^ sin/, sinX,=sin^,sin/, ; 
, . (Tsin I, — (Tsin 1 ■ . . 

hence, smce p= . ' > Pi= t -' we obtam 

^ sin 1 ^ ' sin i 

''■\/»2+»?='W'^j(~2— ^s) (sin^ >t-|- sin^ >.;— 2 sin X sin \ cos %)*. 

This gives the absolute angular velocity with which the principal plane is at any 
instant changing its direction in space ; it is evident that (if the supposition ri=r be 
excluded) it can never vanish except when both planets are in the line of nodes. The 
direction of the rotation is determined by the signs of &>o and »,, 

96. The system of differential equations given in art. 93. affords an example of the 
so-called "elimination of the nodes" in the problem of three bodies. Jacobi, by a 
very remarkable and ingenious transformation, has effected the elimination in a quite 
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different manner*. The equations of art. 93. are merely transformations of the 
original differential equations of the problem, without any integrations ; they are 
however in a form which might perhaps be used advantageously in certain cases for 
the purposes of physical astronomy. Those of Jacobi are obtained by employing all 
the four usual integrals of the problem, and are shown to include an additional inte- 
gration. They have however the disadvantage of substituting the coordinates of two 
fictitious bodies for those of the actual planets, and would probably be inconvenient 
for ordinary practical use ; though in a theoretical point of view they seem to deserve 
more attention than they have hitherto received. It would be wrong to take leave 
of this celebrated problem without referring to another transformation by M. Ber- 
TRAND-^-, which, as has been remarked by a recent writer in the same journal, effects 
six integrations, and therefore represents the furthest advance which has yet been 
made towards a rigorous solution. 

Appendix A. 

When the method described in Theorem VII. (art. 49.) is applied to the solution 
of a system of equations of the form (I.), of which n integrals, a^ ... «„, satisfying 
the conditions [a^, a,]=0, are given, the first step is to express the «-|-l partial dif- 
ferential coefficients ^j &c., and -^; namely, the values of ?/,, ...y^ and — Z in 

terms of «„ ... x„, «„ ... a„ and t. The direct process is then to find X by integrating 
the expression y^dx^-^-y^xA- "-\-yndx^—'Ldt, and afterwards to form the remaining 

integral equations ^=6i, &c. : when this process is adopted, the inferior limits in 

the integrations are perfectly arbitrary ; in other words, we may add to X an arbitrary 
function of a,, a^, ... a„, without altering any of the general properties of the final 
system of integrals. 

But it is generally much more convenient to perform the differentiations with 
respect to a„ ... a^ first, and integrate afterwards; thus we obtain the remaining 
equations in the form 



*.=J(l;''-.+l'''--+"+t'''^--?*)- 



When this plan is followed, the limits are still arbitrary if it be only required that 
the equations thus obtained shall be true; but if it be required that they shall form, 
with the given integrals, a normal solutimi, it is necessary to take the limits in such a 
manner that the functions equated to fi„ \, ... b„ shall be the partial differential coef- 
ficients with respect to a^, a^, ... «„, of one and the same function; which will not 
generally be the case unless care be taken that it should be so. 

In practice, the expression for rfX usually consists of several terms, of which each 

* Comptes Rendus, 1842, part 2. p. 236, &c. t Liouville's Journal, 1852. 

3 A 2 
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contains one of the variables only. Suppose one of these terms is 

p{x, a„ fla, ... a„)dx, 

so that, so far as this term is concerned, we have 

X=4 <p{x, «!, a^, ... a„)dx, 

where A is an arbitrary function of a^, ... a„. Consequently 






and we see that we should not in general obtain the differential coefficients with 

respect to a,, ... a„ of one and the same function X, by merely integrating -7^5 -f-t &c., 

with respect to x, from the same inferior limit A, chosen at hazard. 

But it is evident that we shall attain this end if we adopt the following simple 
rule : — 

Integrate -7^? &c. with respect to x, taking the same inferior limit in each ca§e, 

namely, either 

(1) a value A of x which satisfies the equation (p(x, a,, ... aj=0, or 

(2) any determinate constant (i.e. not a function of a,, ... aj. 

For example, in the problem of central forces (Part I., art. 28, &c.), we had (see art. 19.) 

dX= - hdt-^cdd+(2m{h+(p{r))- k^r'^dr^ {K'-c' sec' xf'dK 
(where r, 6, X are the three variables). 

The very troublesome process of differentiating X with respect to h, k and c after 
first finding X by integrating the above expression, is avoided by the method adopted 
in art. 29 ; namely, by differentiating first, and integrating afterwards. In the inte- 
grations with respect to r, the inferior limit is one of the roots of the equation 

2m{h+p{r)) — Fr-'=0, 

namely (in the case of elliptic motion), the perihelion distance ; and in those with 
respect to X, the inferior limit is ; so that the rule above given is observed. 

At the time of writing the article referred to, neither the rule itself, nor the neces- 
sity of attending to the limits, had occurred to me ; it was therefore, strictly speak- 
ing, accidental that the final integrals were obtained in a normal form. 

In treating the problem of rotation (Section III.), I perceived the necessity of cau- 
tion as to the limits, if the former order of proceeding were adopted ; but preferred 
avoiding the risk of error altogether, by performing the integrations first, so as to 
obtain the actual expression for V. The final equations (R.), art. 44, might however 
be obtained in a more simple way by differentiating_^r*^ ; thus we should have (see 

equations (45.), (46.)), observing that -jj=-j, &c. (art. 44.), and putting 

tt/6 Ai/l 
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(1 — cos^ i~ cos^_y+2 cos i cosj cos 6— cos^ ^)*=Q, 

^.==kU+('^^'^?S^d6]=<. ......... (i.) 

.if. Tcos i cos 3 — cos / 7 .1 « 

In order to get rid of the troublesome integration involved in the term \r%r5«^^j we 

may (1) eliminate this term between the first and last of these equations, and (2) 
eliminate -v^ between the first and second. We thus find the two following equations, 

, Tcos iocs ^— cos; ,. /3 C — A, . /, , v /.. v 

^+i Qsina ^ i^=i-^;C-Acosj.(^+r) ....... (n.) 

.f , /fcosicosfi— cos; ,^1 .fsinflJ 2A,, , , 
which last, combined with the preceding, gives 



MCOSi 



J 



sin W9 k ,, , . a cos « + /3 cos/ .... ^ 



and we may take (i.), (ii.) and (iii.) as expressing the solution of the problem. 

Now we have cos I=?5i!ZL£2?Z|^, ffom which it is easy to find (observing the 

conditions which determine the sign of Q) 

,, cos e cos 8 — cos ;' ,. i . .i , 
—di= Q~~^ — ^^' ^^^ similarly, 

J J cos_; cos S — cos e" 

, ,_, sinWS 

we have moreover a©: 



a 

All the integrations may therefore now be performed immediately; and we may 
take for the inferior limit of 6 any value which satisfies the equation Q=0, or 

(cos ^— cos i cosjy— sin^ i sin^j=0 ; 
this is satisfied by 6=i-\-j, which evidently corresponds to 1=0, J=0, 0=0, and it 
is manifest that equations (i.), (ii.), (iii.) will thus become identical with equations 
(R) of art. 44. 

I do not regret however having introduced the rather prolix investigation of arts. 
39 and 43, because it is interesting to know the actual value of V (equation (48.)), 
which the method just given leaves undetermined. 
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Appendix B. 

On the subject of the transformation of elements, the following additional remarks 
will hardly be superfluous. Suppose Q is oi'iginally a function of the elements a,b, c, 
&c. with t ; and let a, (3, y, &c. be other quantities connected with a, b, c, &c., by 
equations such as 

da=Adci+Bd(5+Cdy+ ...+Kdt, (a.) 

where A, B, C, ... K are given functions of a, jS, y, .... t; or by equations such as 

dtt=Aida+Bidb-\-...-{-Kidt, (a.) 

where A,, &c. are given functions of a, b, ...., t. In either case, if each of the equa- 
tions be integrable per se, we may consider a, b, c, .... as functions of a, (3, y, ...., t; 
and such equations as 

'di=^l^+^Tb + ("•) 

are both significant and true. 

But if the expressions on the right of the equations («.) be not differentials per se, 

the equations (Q.) are either unmeaning or untrue. For the symbol -j- implies one 

of two things; either, that Q is expressed in terms of a, j3, ... ^ without arbitrary 
constants (i. e. that the transformation of Q can be actually effected without inte- 
grating the differential equations of the problem), which is manifestly impossible, 
unless (a.), &c. be integrable per se ; or else, that the differential equations are to be 
conceived to have been completely integrated, so that a, b, &c., and consequently 
Aj, B,, &c., are known as /unctions oft and arbitrary constants, whereby the right- 
hand side of (a.) becomes an explicit function of t (and arbitrary constants), so that 
a, |8, &c. may by integration be expressed in the same way, and, by means of (a.), 
a, b, &c. may be similarly expressed, and finally, by algebraical elimination, a, b, &c. 

become functions of ot, /3, &c., t, and arbitrary constants. On this supposition, ^ 

doc 

has- a meaning, but the equation (Q.) is untrue; for we must have 

da~ da dx'*"db da"^ ■"• ' 
and it is manifestly not true that ^=A, &c. in this case, because the equation 

da= Ada-j-Bd^ -{-... +Kdt, 

not being integrable per se, only subsists for those variations of a, |3, &c. which 
actually take place during the instant dt ; whereas the equation 

da=-^doc+{^dfi+...-\-~dt 

subsists for arbitrary variations of all the variables. This view of the subject entirely 
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coincides, in substance, with that taken by Jacobi ; but the above mode of stating it 
may tend to make it clearer, and to call attention to a matter which, so far as I know, 
is not so much as mentioned in any of the elementary works usually in the hands of 
students of physical astronomy. 

[Jddifion to Appendix B.] 

Received March 1, 1855. 

The remark made above, that the symbol -y- is unmeaning in the case considered, 

is not of course intended to imply that a meaning- may not he given to it ; but then 
such meaning is different from the ordinary signification of the symbol, which is a 
partial derived function. 

The whole matter may be strikingly illustrated by a simple example. 

Consider the movement of a I'igid body about a fixed point. Adopting the nota- 
tion of art. 40 (Part I.), we have 

'pdt-= —cos pd6~ sin ^ sin (pd"^ 

qdt=^ sin ^d6— sin d cos (pd-^ 

rdt=d'p-\- cos dd-4y. 

Let a, |3, y be three new variables, defined by the equations 

nt Pt nt 

a,=z\pdt, ^^=\qdt, ^=1 rdt ; 

Jo Jo Jo 

so that da=pdt, &c. Then the above equations give 

dd= — cos ipdcc-^- sin <pd^, d(p=d'y-{- cot ^(sin (pdcc..-\- cos pd^) 

d-4, = — cosec^(sin^rfa+ cos pd(5). 

Here «, |3, y are the sums of the elementary angles described about the axes in the 
course of the motion ; and no one would maintain that 6, p, 4^ are functions of a, (3, y, 
for the values of the latter variables at any time do not determine the values of the 
former. If therefore we choose to write such equations as 

dS d& . . 

^=-cos<p, ^=smp,kc., 

we must admit that-^j j^, &c. are not partial derived functions in the ordinary sense. 

At most, -J- is the derived function of that function of a which 6 would become if /3 

and y were maintained invariable, i. e. if the motion were restricted to a rotation 

about the A-axis. Again, if we admit such symbols as -r-—-, -— -7r, we must inter- 

° *^ dpdci dadp 

pret them as follows : 

d'^d d dS _ dcos^_ . d(p 
d^—d^ d^— ~ "1^— ^'" ^4' 
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but ^= cot 6 cos (p, and therefore 

jg^= cot ^ sin ^ COS ^ ; 

and in like manner we find the same value for -^-—rj so that in this particular case 

amp 

the condition ^=||^ is verified. 

AH d^S 

But if we take j-^ and j-j- in the same way, we find the former = sin (p and the 

latter =0, so that in this case the condition is not verified. The geometrical meaning 
of this is obvious ; analytically it is merely an instance of a general fact, pointed out 
by Jacobi ; namely, that the effect of two successive pseudo-differentiations with 
respect 'to two independent variables, is not generally independent of the order of 
operations. 

If V be the potential of another body, given in position, upon the body considered, 
then V is a function of 6, p, ip, and 

dy=^dd+^d^+j^d^; 

and if we substitute for d6, d<p, dip their values in terms of da, dj3, dy, we obtain an 
expression which we may call Lc?a+Mc?j3+N6?y, L, M, N being functions of ^, (p, -p, 

-17' -T-» -TT' of which, as is well known, the mechanical meanings are the moments of 

d$ d(p d^ ' ' '^ 

the attraction of the second body about the three axes. Here again no one would 

dY 
maintain that V is a function of a, /3, y ; and if, as is often done, we say -j-=L, &c., 

the above remarks apply in all respects to these equations. 

I should have thought it superfluous to dwell so much on these points if it had not 
appeared that writers on physical astronomy have in some instances either overlooked 
the distinction between real and pseudo-differentiation, or at least have failed to point 
it out to their readers. The only discussion of the subject which I have met with is 
that given by Jacobi, in the correspondence referred to. 

7TT 

It may be added, that in general investigations, where symbols such as -j-^ &c. 

may be used without defining the nature of V, or the precise meaning of a, /3, &c., 

dW d^Y 
serious errors might be committed if it were assumed that the condition -7—j:;=-5s-t- 

° map dpaa 

always subsisted. 



Appendix C. 

The theorems relating to the transformation of coordinates, given in Section VI., 
may be made more general, and in many cases more useful, as follows : — 
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If Xi, j?2, ... x^ be the coordinates employed in the first statement of any dynamical 
problem, the differential equations are comprehended in the formula 

^.{(f)'-^^^-O (D.) 

[If there be any forces, such as those aiising from a resisting medium, which do 
not satisfy the natural conditions of integrability, then on the right-hand side of the 
formula (D.), instead of we shall have an expression such as ^^(XjS^i) ; but such 
terms are easily introduced and allowed for separately, and do not affect the follow- 
ing investigation. I shall therefore here assume that they do not exist.] 

In the above formula, W is a given function of x^, ... x^, x\, ... x'„ which may also 
explicitly contain t. 

In Section VI. the only case contemplated was that in which x^, ... x„ are inde- 
pendent coordinates ; in which case tiie formula (D.) is equivalent to n separate equa- 
tions, since Ix^, &c. are wholly arbitrary and independent. 

In practice, however, it is often more convenient to use, at first, a set of coordinates 
more in number than the independent variables of the problem, and therefore subject 
to equations of condition. 

Let us assume then that the n coordinates x-^, ... x„ ate suiyect to r equations of 

condition, 

L,=0, L2=0, ... L,=0, 

where Li, &c. may explicitly contain t, besides the n variables ^i, &c. 

If we introduce the n conjugate variables y^, ... y^ defined by the equations y(=-T-T' 

i 

and take Z a function of «„ &c., ?/„ &c. (with or without t), defined by the equation 

z=-[w]+[x;>,+...+My» 

(the brackets indicating thato^'i, &c. are expressed in terms ofy„ &c.), then it follows 
exactly as in art. 18 (Part I.), that the formula (D.) will be changed into the system 

dZ 



Xi=z- 
dyi 



K^^+S)^^^=o. 



(E.) 



[In the most usual problems W is of the form T+U, where T is homogeneous and 
of the second degree in x\, &c., and U does not contain x\, &c. at all. In this case 
Z is only T— U expressed in terms oly^, &c., instead of x\, &c. But T is not neces- 
sarily homogeneous ; in fact it is not so in problems relating to motion relative to 
the earth, as affected by the earth's rotation,] 

Let us now suppose that the system (E.) is to be transformed by the introduction 

of the m independent coordinates li, fj, ... |„, and of the new conjugate variables 

ni, I2, ••• fimi where m=n—r. And let it be required to investigate a theorem by 

means of which the transformation may be effected without recurring to the miginal 

formula (D.). 

MDCCCLV. 3 B 
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The definitions of the new coordinates |i, &c. will furnish m equations (which may 
explicitly contain t) by means of which |„ ... |„ may be expressed as functions of 
Xj, ... x^ (with or without t) ; and conversely, by means of these m equations, together 
with the r equations of condition Li=0, &c., the n variables x^, x^, ... x„ may be 
expressed as functions of fi, ... |„, with or without t. When Xj, ... x^ are so expressed, 
let them be represented by [Xi), ... {x„). We shall have then 

so that Xi, &c. are expressible (and in only one way) as functions of Si, &c., Ii, &c. 

If then the formula (D.) be transformed by expressing x„ &c., x\, &c. in this manner, 
it becomes, as is well known, 

where (W) represents the result of transforming W as above ; and since S|i, &c. are 
now independent, this formula breaks up into the m separate equations 

diV\[) 
Moreover, if we now define % by the equation -W^=?;i, and put 

T=_(vv)+(r.K+...+(CK, 

where (li), &c. are expressed in terras of yii, &c., we know already (art. 18.) that the 

system (F.) becomes 

d^ ,_ d^ 

^*=^' ''^--^ ^^-^ 

Now let P be a function of the m new variables li, ... |„, and of the n old variables 
J/i! ••• ^n (with or without t), defined by the equation 

P= (^i)yi + (■^2)^2+ — + (««)«/«. 

Since ;jj=:-W- , and since (W) contains |-, &c., only through x-, &c., we have, ob- 

, dW 
serving that -^=3/;, 

but from the equation {x\) we have ip-—m"' 

consequently rti-y^-^-iry^-^ -\-...-{-y^-^, 

an expression evidently equivalent to -^' Thus ^j, may he defined hy the equation 

dV ^ 
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without recurring to the formula (D.). And since each equation of condition gives, 

dZ 
dyi 



All 
if differentiated totally with I'espect to t, with the substitution of t- for «■, an equa 



dt 



tion such as 

^— dt'^dw^dy^'^'-'^dx^dyr, ^■^•) 

the m equations (;?.) with the r equations L (in which last x^, &c. must be expressed, 
after the differentiation, in terras of |, &c.), give n equations by means of which 
r/i, ...yn may be expressed in terms of li, &c,, n-^, &c., and can be so expressed only in 
one way. 

Lastly, the value of "^ (see equations (G.)), may be obtained as follows: — 

Since •^F=_(W) + (i;);;,+ ... + (Q^„ 

and Z=-[W] + [^;]3/:+... + [4]y«, 

and (W) is only [W] differently expressed, we have, without reference to modes of 

expression, 

^-Z=2,(&)-E,(%). 

On the other hand, since P can contain t explicitly only through {x^, &c., we have 

'_ / d{w>) \ 

'—^^\yi dt )> 

and also 2,(%)=2,(z/,^) +S,(f g) ; 

dP 

hence, observing that ^.=>ii, we obtain 

, , dV 

dP 

and therefore, finally, Y=Z— -^' 

dP 
so that Z—-^ must become identical with W, when expressed entirely in terms of 

the new variables. 

These results may be stated in the form of the following fhecyrem. The system (E.) 
is transformed into the system (G.) by the following substitutions : — 

(1) a?„ ... x^ are expressed in terms of |„ ... |„ by means of the m equations which 
define the latter variables, together with the n—m equations of condition 

Li=0, ... Ly=0 (r=n—?n). 

dP 

(2) ill, ... ;?„ are defined by the m equations ■Ty-=ni-, where the modulus of trans- 
formation P is given by the equation 

P= (*Oyi + (^2)^2+ . . •+ (^»)3/«, 
(*i), &c. being expressed in terms of li, &c., so that P is explicitly a function of g„ ... |„, 
y, ...yn, with or without t. 

2b2 
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dF 
(3) If is defined by the equation 'V=Z — -jT} in wliich (after the explicit differentia- 
tion of P with respect to t), x,, &c., i/,, &c. are to be expressed in terms of the new 

dP 
variables, y^, &c. are thus expressible by the help of the m equations -^=;;jand the 

dh . „ rdJ^d7i\_ 



du , „ /dL, d'L\ 
n-m equations ^+2,(^ ^.j=0 

If {xi), &c., do not contain t 
expressing Z in terms of the new variables 



dV 
If {xi), &c., do not contain t explicitly, then -^=0, and Y is obtained merely by 



It may be observed that the whole of the above reasoning would apply to the case 

in which the new variables |i, ... 1^ are more in number than the independent vari- 

dP 
ables of the problem (or m>n—r), with this exception ; that the m equations jp-=7j;, 

together with the r equations obtained by differentiating the equations of condition 
totally with respect to t, would be more than sufficient to express ?/i, ... y^ in terms of 
the new variables ; consequently ?/„ &c. might be so expressed in different ways, and 
therefore, although the value of ""F obtained by the above rule would certainly be the 
same as that obtained by recurring to the original formula (D.), the. form of Y might 
be different, and therefore the resulting formula erroneous. 

There must doubtless exist some rule for choosing n — m combinations of the 
equations of condition in such a way as to lead to the correct forms of y^, ... y„ as 
functions of the new variables ; but I have not at present attempted to invQ^stigate it, 
and perhaps it would be hardly worth while. The theorem in the case in which the 
new coordinates are independent, may, I believe, be practically useful. 



ERRATA IN PART I. 

Art. 1. equation (4.), for dx read dxi. 

Art. 10. In paragraph preceding equation (26.) omit the words " not containing t explicitly." 

Art. 18. equation (/3), for yi read y^. 

Art. 19. equation (29.), for h read hi. 

Art. 24. second line after equation (L.),/o»' " such as h, k" read " such as/, g." 

Art. 30. The expressions equated to h, k, c, and the three terms in the left-hand column of the table of 

elements, should each be multiplied by m. 
Art. 42. near the end, for " according as 9 is between o and ir, or not" read " according as is between 

IT and '2ir, or between o and tt." 



